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Institute of Physics, Polish Academy of Sciences, Warsaw 


(Received June 6, 1960; translated paper received on August 6, 1960) 


The formation of etching figures on the (0001) face of bismuth selenide single 
crystals is investigated. Different etching mixtures were found to give rise to figures 
of different shape and density of distribution. The mixture best adapted for the detection 
of dislocations was that (No. 5) containing bromine and iodine. In order to explain the 
different forms and densities of the figures, depending on the procedure of etching, a mod- 
el is proposed wherein the walls of the etching figures are considered to consist of an 
assemblage of steps. The model relates the thickness of the steps to the form and orienta- 
tion of the figures and to the rate of etching. 


Recently, a number of papers on methods of revealing details of the defect 
structure of crystals have been published. 

One of the most involved and least elucidated problems is that of dissolving the 
crystal. This process often leads to widely different results according to the defects 
present in the crystal, the type of structure and the method used. Maybe the etching 
figures formed at the intersections of dislocations with the crystal surface are the 
most readily interpreted. 

The present investigation aimed at interpreting problems of this type in the case 
of semiconducting bismuth selenide crystals of purity corresponding to 8 x 1018 — 
— 2 x 10% holes per cm? (type p), etched chemically. 

In order to obtain insight into the process of formation of various figures according 
to the conditions of dissolving, a number of etching substances were utilized. Thus, 
the effect of small admixtures to the etching solution on the form of the figures ob- 
tained on the crystal surfaces was detected. 

The experimental procedure is the following: Plates of bismuth selenide obtained 
by splitting along the (0001) plane were immersed in the etching solution for two 
minutes, rinsed in ethanol and dried. The plates were etched at room temperature; 
if etching was continued, the temperature of the crystals rose and the figures obtained 
in the process of etching were not easily reproducible. The preparates thus obtained 


(103) 


104 


Picture 2 


were subjected to observation with a metallographic microscope and micro-interfero- 
meter magnifying 120 to 1400 times. 
The form of the etching figures varied according to the solution used. The com- 
position of the latter is given at the end of the present paper, for the various cases. 
If mixture No. 1 was used, the figures obtained did not exhibit the symmetry 
characteristic of the erystal face considered. 
With mixture No. 2, a more distinct picture of the figures was obtained; this 
consisted of triangular pits whose. edge was about 7 long, and which were distributed 
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Picture 3 


Picture 4 


with a density of the order of 108 per cm?. The pits were flat and, in some regions, 
were arrayed in rows. The formation of the pits proceeded as long as etching went on, 
at first presenting a triangular shape and then rounding, finally to become diffluent. 
The density of the pits having sharply defined edges was constant throughout etching. 
This would seem to indicate that these figures derive from defects of rather small 
linear dimensions and lattice irregularities distributed evenly throughout the crystal. 

Picture 1 presents single etching figures of different shapes possessing or not the 


triangular symmetry. 
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Picture 5 


Picture 6 


Picture 2 shows etched slip bands, this fact was established as the result of 
complementary X-ray studies. 

Some regions exhibit rows of etch pits forming a lattice of triangular lattice 
points (Picture 3). Most probably, these are boundaries between blocks rotated mu- 
tually by a small angle. 

Etching mixture No. 3 yielded faces covered with triangular figures outlined 


sharply. Picture 4 shows the (0001 crystal face on etching with this mixture. Some 
p^ à g 
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Picture 7b 


figures are seen to be arranged in larger ones of the same shape and symmetry. Picture 
5 shows etch pits grouped in rows. 

Mixture No. 4 yielded etching figures resembling though somewhat more contrasty 
than those obtained with No. 3: the single figures were very sharply outlined, whereas 
the groups of small figures covered larger areas and partook of the character of large 
cavities of complex structure reminiscent of the figures observed by Oberly on german- 
ium (Oberly 1957). É 

Picture 6 is that of the surface of a Bi Ses crystal subjected to etching with mixture 


No 4, showing either kind of figures, 
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Pictures 7a and 7b show the structure of the large cavities of Picture 6 at greater 
enlargement. 

Mixture No 5. also yielded two kinds of figures having sharp contours and the 
symmetry corresponding to the erystal face subjected to etching: the large ones 
appearing sporadically, and the small ones distributed densely. The large pits resulted 
from deformation in the crystals investigated; they sometimes appeared in rows 
(Picture 14). 

X-ray investigation according to Dr J. Auleytner’s method (Auleytner 1958) 
and utilizing a point source of X-rays with focus diameter of about 20u (Au- 
leytner 1957) established that the large pits appear at the boundaries of mutually 


Picture 8 


rotated regions, corresponding to sites at which the dislocations intersect with the 
crystal surface. 

Picture 8 shows a crystal surface covered with a group of dislocation etch pits, 
whereas Picture 9 is that of a dislocation boundary terminating in a group of disordered 
dislocations. 

In order to investigate more closely the surface etched in mixture No. 5, a number 
of photographs were taken utilizing a microinterferometer. This led to a more exact 
determination of shape and structure as well as of the depth and angle of inclination 
of the walls bounding the pits. The walls of the pits formed an angle ranging from 
6° to 9° with the crystal surface, the maximum depth attaining about 6u, whereas 
that of the small pits amounted to about 0.4 It. 

Moreover, the pictures obtained with the microinterferometer served to establish 
the depth of the inequalities on the surfaces subjected to etching, in the regions present- 
ing no pits of either kind. Fluctuations of the inequalities of the surface between 
cavities amounted to some tenth of a micron, 


109 


Picture 9 


Picture 10 


Picture 10 is that of a region covered with small etch pits. 

Picture 11 shows separate large etch pits arising from dislocations presenting 
different angles of inclination with respect to the crystal surface of the walls bounding 
them. Hence it is to be inferred that the dislocation lines giving rise to the pits were 
inclined with respect to the crystal surface at an angle other than 90°. 

When observing the distribution of the interference fringes on the surface of the 
walls bounding the etch pits, it is seen that the angle at which they are inclined with 
respect to the crystal surface varies. Hence it is to be concluded that the dislocations 
giving rise to the cavities in the process of etching were of varying inclination 
throughout a layer whose thickness equals that of the cavity. 
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Picture 12 


Picture 12 shows large etch pits of flat bottom, deriving from dislocations. This 
may result from abrupt curvature of the dislocation lines (Yamamoto 1958). 

Picture 13 shows some groupings of pits corresponding to dislocation lines. 

Picture 14 represents an etched dislocation boundary. The regular distribution of 
the dislocation pits along the boundary is clearly visible. The interference fringes 
point to an angle of inclination of the latter with respect to the crystal surface differing 
from 90°. 

In all cases, the surfaces bounding the pits were found to belong to the same 
type of crystallographic form. 

In considering the process of etching, the rate of etching is sometimes represented 
as the sum of a component. normal and of one tangent to the surface subjected to 
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Picture 14 


etching (cf., e.g., Yamamoto 1958, Gilman et al. 1958). This procedure is especially 
well adapted to the case of the (0001) face of bismuth selenide crystals, as the latter 
present a layer structure involving considerable anisotropy of their properties. This is 
especially salient when the c-axis direction is compared with the directions perpendic- 
ular thereto. 

The rate of etching depends on the energy of the bonds retaining the atoms on 
the crystal surface. If a step is formed on the (0001) face, the rate of etching in the 
tangential directions will greatly exceed that in the normal direction, since the atoms 
on the edge of the step are much more loosely attached to the crystal than. those on 


the (0001) face, 
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Hence, the etching rate may be said quite generally to depend on two factors, 
namely: 

1) the rate of formation of monatomic pits, and 

2) the etching rate in the tangential directions. 

Either factor depends on the energy per single bond of the atoms forming the 
crystal surface. The reaction rate of etching is, among others, a function of exp (E/kT), 
wherein E denotes the activation energy dependent on the energy of the bonds of a given 
lattice point on the crystal surface. 

Taking into account the ratios of the etching rates for the various possible crystal 
faces, it is obvious that, as the temperature rises, the rates of etching of the different 
faces and edges will tend to assume one value. Thus, a rise in the activation energy 
will, on the one hand, diminish the etching rate and, on the other, enhances the 
anisotropy of etching. 

Detailed consideration of the structure of the Bi,Se, crystals investigated yields 
an explanation of the different forms and densitites of the etching figures depending 
on the method of etching employed. Bismuth selenide crystallizes in the trigonal 


system, in the space group D$, = R3m (Semilyetov 1954). The basis of the Bi,Se, 
crystal structure is given by nine layers of selenium atoms disposed according to 
the symmetry of the trigonal system 


r.ABAGBGAGA.- 
wherein the bismuth atoms occupy two thirds of the octahedric interstices 
AcBAcBAcBaCbACbA... 


Taking into account the layers transforming into themselves by translation, we have 


ABCA ... 


On the (0001) face of erystals possessing such structure, etching can lead to the forma- 
tion of two kinds of differently oriented pits. The pits of either kind have edges wherein 
the density of atom setting is the greatest. 

In order to distinguish pits of different kinds and the same configuration, one 
kind will be denoted by Z, and the other by the letter P. 

Fig. 1 represents a single layer, showing the direction of the bonds with the next 
and the preceding layer wherein two pits are present: one of configuration Z, and the 
other one of configuration P. In the circles denoting atoms or molecules, blank fields 
represent bondings with the upper layer, black ones — with the lower layer, and 
bondings with neighbouring atoms within the layer are.symbolized by the boundaries 
between the black and blank fields. From Fig. 1, pits of either configuration are seen 
to be equally possible if they are of monomolecular depth, as then the atoms on their 
edges are equally strongly bonded to the base. If, however, the etch pits become 
deeper conditions can arise when their walls will consist of steps many layers thick. 

According to the structure of the pits’ walls, the etching rates for cavities of 
either configuration can undergo differentiation or can remain equal. This depends 
on the distances and thickness of the steps bounding them, 

The following cases can occur: 
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1) the pits can be bounded by mutually separate steps of monomolecular thickness, 

2) the pits can be bounded by mutually separate steps of thickness exceeding 
that of one layer (Fig. 2), or 

3) the pits can be bounded by surfaces exhibiting no steps. 

In the case of (1), as in that of the monomolecular pit, both configurations are 

- equally possible. As a results of this, the etching figures can assume a macroscopically 
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Fig. 1. Two equally probable antisymmetric etching figures on the (0001) face, of depth equalling the 
thickness of a single layer. 

The circles represent atoms: blank fields denote bonds directed upwards, black fields denote bondings 
with the neighbouring lower layer (omitted for clarity), and the black/blank boundaries denote bondings 
with neighbouring atoms in the same layer. 

The atoms on the edges in either figure (L and P) are seen to be equally strongly attached to the crystal. 


more or less regular shape not related to the crystallographic symmetry of the surface 
wherein they are present. This case can serve for describing the figures obtained in 
etching with mixture No. 1. 

In the case of (2), the pits obtained by etching would present walls constructed 
of mutually separate steps of thickness greater than that of a single layer. Fig. 2 shows 
two pits of configuration L and P, of the depth of one two-layer steps. 

The bonds attaching the atoms of the edges of the upper layer forming the upper 
surface of the step to the substratum are seen to be equally strong in the case of either 
pit; however, the atoms on the edge of the next layer in either pit are bonded to the 
substratum with different strength (it should be held in mind that, on the circles 
denoting atoms or molecules, blank fields represent bondings with the upper layer, 
black ones — with the lower layer, and bondings with neighbours in the same layer 
* are symbolized by the boundaries between black and blank fields). 
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From Fig. 2, the atoms on the edges of the layer immediately beneath the upper 
one are less strongly attached to the substratum in the case of configuration Z than in 
that of configuration P. In case L, the edge of the layer underlying the upper one 
contains atoms having two non-saturated bondings with neighbouring atoms and 
two with atoms of the upper layer, whereas in case P the same atoms present two 
non-saturated bondings with neighbouring atoms and one with those of the upper 


layer. Hence, the etching rate for the lateral walls of the steps in configuration Z should 
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Fig. 2. a) Pits of the same configurations as in Fig. 1, albeit of increased depth and exhibiting the forma- 
tion of a step of the thickness of two layers. 
The notations is that of Fig. 1, the atoms of the upper layer (forming the crystal face etched) being re- 
presented by blank fields (I), the atoms of the following layer (II) — by blank shaded fields, and the blank 
fields of the atoms in the layer constituting the bottom of the pit (III) are dotted. 
It is seen that, although the atoms on the edge of the upper layer are equally strongly bonded to the 
substratum in either configuration of the pit, those on the edge of the underlying layer are less strongly 
attached in configuration L than in configuration P. This explains why the crystal, on etching, exhibits 
pits of configuration P only. 
b) Section of L cavity. 
c) Section of P cavity. 


exceed that for configuration P, and the etch pits should be expected to assume con- 
figuration P. Moreover, it is obvious that, as the thickness of the steps increases, so 
will also the differences between the etching rates for the lateral walls of the steps in 
either configuration of the pits. Simultaneously, the etching rate of the lateral walls 
of the steps in either configuration will diminish, since then the density of step edges 
per unit depth of the figure etched diminishes, and the atoms on the edges are those 
that are the least strongly attached to the substratum. 

Case (3) can be said to represent the limit for case (2) with the pit walls constituting 
a homogeneous step surface and the three walls being crystallographic planes inclined 
at a constant angle with respect to the crystal surface subjected to etching. 

From the model proposed it may be concluded that, as the etch pits assume forms 
that are progressively sharper, the etching rate of the crystal should decrease, and, 
moreover, that all cavities should have the same orientation, as was indeed the fact 
when applying mixtures Nos. 1—5. 

At the same time, further observation established that the pits should possess 
a structure consisting of steps of a thickness of several layers, and not be bounded 
by surfaces devoid of steps. This results from the interferograms, which point to pits 
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bounded by surfaces of variable inclination with respect to the erystal surface, proving 
that case (3) cannot occur. Obviously, the model proposed, wherein the wall is consider- 
ed to consist of an array of steps, can also be adapted to crystals belonging to systems 
other than bismuth selenide. 

From the variable inclination of the surfaces bounding the etch pits, the dislocations 
pervading the erystal can. be concluded to have been curvilinear (Yamamoto 1958). 
However, the angle of inclination of the dislocations accounting for the formation 
of the pits cannot be determined unless the angular distribution of the etching rate 
is known. 

Investigation of the effect of admixtures of bromine and iodine to the etching 
liquid suggests that bromine and iodine raise the activation energy, thus leading to 
a differentiation of the size of the etch pits, and that they may enhance the formation 
of multi-layer steps. 


Results 


1) The effect of various etching mixtures on the microscopic structure of the etched 
surfaces of bismuth selenide single crystals was investigated. 

2) Etching mixtures Nos. 1—5 were found to produce different etching figures, the 
density and shape thereof varying according to the mixture used. 

3) The addition of bromine and iodine to the etching mixtures was found to produce 
etching figures possessing sharper outlines; moreover, the figures, then fall in two 
classes, namely, densely disposed minute pits, and large ones of a dislocational 
origin. Hence, bromine and iodine can be said to be factors influencing the contrast- 
iness of the picture. 

4) Mixture No. 5 was found to yield the best results in revealing dislocations. 

9) An explanation is proposed for the fact that a single type of pits is observed on 
the crystals etched, on the basis of the structural data for bismuth selenide. 

6) The mechanism of etching is discussed on the basis of the experimental results 
obtained and with respect to the model proposed, which is applicable to crystals 
of other systems also. 


Etching mixtures applied: 


No. 1: 2 parts of HNO, No. 2: 2 parts of HNO, 
l part of HCl l part of HCl 
l part of H,O, l part of H,O, 
No. 3: 2 parts of HNO, No. 4: 2 parts of HNO, 
1 part of HCI l part of HCI 
Br, to saturation l part of C H OH 


Ją to saturation 
No. 5: 2 parts of HNO, 


l part of HCl 
1 part of C,H,OH 


Br, and J, to saturation _ 
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EINFLUSS VON FLÜSSIGKEITEN AUF DIE TRIBOERREGUNG 
DES ALUMINIUMS ZUR PHOTOSTIMULIERTEN 
EXOELEKTRONENEMISSION 


Von T. LEWOWSKI UND B. SUJAK* 
Institut für Experimentalphysik der Universität, Wroclaw 


( Eingegangen am 7 Juni, 1960) 


Der Einfluss von Flüssigkeiten auf die Triboerregung des Aluminiums zur Exoelektro- 
nenemission wurde untersucht. Die Intensitit dieser Emission wurde bei Beleuchtung 
mit weissen Licht einer Wolframlampe mit Hilfe eines offenen Luftspintzenzählers ge- 
messen. Die untersuchten (vorwiegend organischen) Flüssigkeiten zeigten einen unter- 
schiedlichen Einfluss auf die Triboerregung je nach der Erregungsart: Bürsten mit einer 
Stahlbürste oder Abschabung mit einer Feile (auch Messerkante). Im Falle der Bearbeitung 
der Oberfläche mit einer Stahlbürste lassen sich die untersuchten Flüssigkeiten in drei 
Gruppen einteilen: 1) diejenigen die stürker lóschend auf die Triboerregung wirken als 
die Luftatmospháre, 2) diejenigen die milder lóschend wirken und 3) diejenigen die den 
gleichen Einfluss auf die Triboerregung wie die Luftatmosphäre zu zeigen scheinen. 
Die Moleküle der Flüssigkeiten der ersten Gruppe enthalten sowohl Sauerstoff- als auch 
Wasserstoffatome, die der zweiten Gruppe nur Wasserstoffatome. In den Molekülen 
der Flüssigkeiten der dritten Gruppe sind weder Wasserstoff- noch Sauerstoffatome vor- 
handen. Im Falle der Abschabung der Oberflüche mit einer Feile oder einer Messerkante 
wirkten dagegen alle untersuchten Flüssigkeiten stets mehr oder weniger stärker lóschend 
als die Luftatmosphäre. 

Um den Einfluss der Flüssigkeiten auf die Triboerregung gut verfolgen zu können, 
musste der Kontakt Aluminium-Flüssigkeit in einer möglichst kurzen Zeit nach der 
Oberflächenbearbeitung stattfinden. Der Einfluss von elektrischen Strömen in den Flüssig- 
keiten auf die Triboerregung der Aluminiumelektroden wurde orientierungsweise auch 
untersucht. Es zeigte sich, dass je nach der gewählten Stromrichtung eine elektrolytische 
Erregung oder Auslöschung der photostimulierten Exoelektronenemission stattfinden 
kann. 

Die chemische Wirkung der an der Kontaktoberfläche sich zersetzenden Moleküle 
wird als Ursache des Einflusses der Flüssigkeiten auf die Triboerregung zur photosti- 


mulierten Exoelektronenemission gedeutet. 


Einführung 


In unserer Arbeitsgruppe wurden verschiedene Oxydationsprozesse an Alu- 
miniumoberflächen mit Hilfe der Exoelektronenemission untersucht. Es wurde dabei 
unter anderen festgestellt, dass durch das Spülen der Proben mit Wasser die Tri- 


* Zur Zeit in The Clarendon Laboratory, Oxford 
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boerregung zur Exoelektronenemission in kürzester Zeit beträchtlich erniedrigt 
werden kann. Diese löschende Wirkung des Wassers, die Benutzung von Ätzverfahren 
bei den Untersuchungen von Exoelektronenemission, sowie das Problem der Elektronen- 
emission in Flüssigkeiten im allgemeinen, haben uns veranlasst die Einwirkung von 
Flüssigkeiten auf die mechanische Erregung zur Exoelektronenemission näher zu 
verfolgen. Die vorliegende Arbeit enthält einen Bericht über die Ergebnisse dieser 
Untersuchung an Oberflächen des technischen Aluminiums in einer schematischen 
Darstellung. Die Ergebnisse scheinen auch für Korosionsfragen der deformierten 
Metalle von Bedeutung zu sein (Thomson 1958). 


Untersuchungsmethode 


Zur Untersuchung der Einwirkung von Flüssigkeiten auf die Triboerregung zur 
photostimulierten Exoelektronenemission (die weiterhin als Ph. Ex. Em. bezeichnet 
wird) wurden Proben aus technischem Aluminium vorbereitet. Die Länge der Proben 
betrug etwa 8 cm, die Breite etwa 2 cm die Dicke 0,2 cm. 

Ein Teil der Proben wurde mit einer Stahlbürste oberflächlich bearbeitet. Der 
andere Teil dagegen wurde mit Hilfe einer groben Stahlfeile oder eines Stahlmessers 
abgeschabt. Die Bearbeitung der Oberfläche erfolgte immer in atmosphärischer 
Luft unter Tageslicht. Kurz nach erfolgter Bearbeitung der Oberfläche wurde die 
Verteilung der Emissionsitensität der Ph. Ex. Em. entlang der ganzen Probe gemessen. 
Dieses Abtasten der Proben geschah in Dunkelheit mit Hilfe eines offenen Luft- 
spitzenzählers durch welchen hindurch das weisse Licht einer Wolframlampe die 
untersuchte Stelle der Proben beleuchtete (Sujak 1957, Mader und Sujak, 1960). 
Somit wurden für jede Probe die Abtastkurven N/t = f(x) der Ph. Ex. Em. aufge- 
nommen. Die benutzte Messanordnung ist in Abb. 1. schematisch wiedergegeben. 

Zur weiteren Untersuchung gelangten nur Proben deren Abtastkurven einen 
gleichmässigen und horizontalen Verlauf zeigten. Diese wurden danach in der ge- 
wählten Flüssigkeit (Raumtemperatur) teilweise für etwa 5 Sekunden eingetaucht, in 


Abb. 1. Schematische Darstellung der Messanordnung: (1) Spitzenzähler, (2) Al-Probe, (3) Bewegliche 
Probenunterlage, (4) Papierband-Kuppelung, (5) Verstärker, (6) Impulsmittelmesser, (7) Hochspannungs- 
speisegerät, (8) Selbstschreibegerät zur Aufzeichnung der Ph. Ex. Em. Abtastkurven, (9) Lichtquelle 
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kaltem Luftstrom (Raumtemperatur) getrocknet und die Abtastkurven der Ph. Ex. Em. 
wurden dann wieder aufgenommen. In sämtliche Abbildungen wurden die Abtast- 
kurven der Ph. Ex. Em., die nach der mechanischen Bearbeitung aufgenommen 
wurden, gestrichelt eingezeichnet. Voll ausgezogen wurden dagegen die Abtastkurven, 
die nach der Einwirkung der angegebenen Flüssigkeit gemessen wurden. 

Um den Einfluss des Abklingens der Emissionsfähigkeit (der Erregung) des 
nicht eingetauchten Teiles der Probe in der Luftatmosphäre selbst oder in den Flüssig- 
keitsdämpfen berücksichtigen zu können, wurde immer eine Kontrollprobe benutzt, 
die in die Flüssigkeit überhaupt nicht eingetaucht wurde. 

Zur Untersuchung kamen verschiedene, grösstenteils organische Flüssigkeiten, 
in deren Molekülen Wasserstoff- und Sauerstoffatome, oder nur Wasserstoffatome, 
oder weder Wasserstoff- noch Sauerstoffatome eingebaut sind. Der Reinheitsgrad der 
Flüssigkeiten war „chemisch rein“ oder „Pro analisi“. 


Versuchsergebnisse 


l. Mechanische Bearbeitung und danach folgendes Eintauchen 
in eine Flüssigkeit 

Die Flüssigkeiten zeigten verschiedene Einwirkung auf die Abtastkurven der 
Ph. Ex. Em., je nachdem, ob die Al-Proben vorher mit einer Stahlbürste oberflächlich 
gerieben oder mit Stahl (Feile, Messerkante) abgeschabt wurden. 

a) Mit Stahlbürste bearbeitete Oberfläche. In diesem Falle können alle unter- 
suchten Flüssigkeiten in folgende drei Gruppen eingeteilt werden. 

1. Gruppe. Zu ihr gehören alle diejenigen Flüssigkeiten, die eine Abnahme 
der Erregung zur Ph. Ex. Em. im Vergleich zum nichteingetauchten Teil der 
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Abb. 2. Ph. Ex. Em. Abtastkurven (schematisch) für Al-Proben, die vorher mit Stahlbürste bearbeitet 
wurden (— — — —) vor dem Eintauchen, (— — — — —) nach dem Eintauchen in eine der Flüssigkeiten 


a) Nach der Einwirkung einer Flüssigkeit der 1 Gruppe (H- sowie O-Atome im Molekül). : 

b) Nach der Einwirkung einer Flüssikgeit der 2. Gruppe (H- aber keine O-Atome E MEREUN. 

c) Nach der Einwirkung einer F lüssigkeit der 3. Gruppe (keine H- sowie keine De ome im Bol 

Die Lage der Al-Probe ist über den Abtastkurven angedeutet. Der mit der Flüssigkeit in Berührung 
kommende Teil ist gestrichelt gezeichnet, 
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Al-Probe bewirken, wie Wasser (H,O), Athylalkohol (C,H,OH), Athylather 
(C,H,OC,H;), Essigsäure (CHCOOH) und Azeton (CH,COCH,)*. Siehe Abb. 2a. 

2. Gruppe. Zu dieser werden die Flüssigkeiten gerechnet, die den zeitlichen 
Abfallprozess der Erregung zur Ph. Ex. Em. verlangsamen. Dies wird als eine Zunahme 
der Erregung zur Ph. Ex. Em. des eingetauchten Teils der Probe im Vergleich zum 
nichteingetauchten Teil gemessen (siehe Abb. 2b). Hierzu gehören: Benzen (CgHs), 
Toluol (CgH;CH,), Chloroform (CHCL,), Trichloráthylen (C HCl). 

3. Gruppe. Zu dieser werden diejenige Flüssigkeiten zugeteilt, die auf die 
Fähigkeit zur Ph. Ex. Em. keinen Einfluss zu haben scheinen, z. B. Chlorkohlenstoff 
(CCH) Siehe Abb. 2c. 

b) Mit Stahl abgeschabte Oberfläche. Für alle untersuchten Flüssigkeiten zeigt 
sich, dass in dem Eintauchsbereich stets eine Abnahme der Intensität der Pb. Ex. Em. 
im Vergleich zum nichteingetauchten Teil eintritt, wenn die Oberfläche der Al-Probe 
mit einer Feile oder Messerkante vorher abgeschabt wurde (Abb. 3). 

Um festzustellen ob bei den beobachten Erscheinungen die elektrische Leitfähig- 
keit der Flüssigkeiten keine ausschlaggebende Rolle spielt, wurde die Leitfähigkeit 
des Azetons mit Hilfe von gelösten KJ um Grössenordnungen geändert. Es konnte 
jedoch kein deutlicher Einfluss der Leitfähigkeit des Azetons auf seine Löschwirkung 
gefunden werden, wie dieses nach den Erwägungen von Schaafs (1959) zu erwarten 
wäre. Wie bekannt, betrachtete Schaafs lokaliesierte, hohe Spannungsgradiente als 
Ursachen der Emission von Exoelektronen nach mechanischer Erregung. 


SS OT TT aT TT 


Abb. 3. Ph. Ex. Em. Abtastkurven (schematisch) fiir Al-Proben, die vorher mit einer Feile oder Messer- 
kante abgeschabt wurden: (— — — —) vor dem Eintauchen, (~~) nach dem Eintauchen 
in eine Fliissigkeit irgendeiner der drei Gruppen. 


2. Mechanische Bearbeitung, Eintauchen in Flüssigkeiten und 
danach folgende Erwärmung. 

Die Al-Proben wurden wiederum zuerst entweder mit einer Stahlbürste bearbeitet 
oder mit einer Feile (Messerkante) abgeschabt. Die Proben wurden danach 5 Sekunden 
lang in einer Flüssigkeit der 2. Gruppe zur Hälfte eingetaucht gehalten, im Luftstrom 
kalt getrocknet und auf Ph. Ex. Em. Abtastkurven untersucht. Die Proben wurden 
danach eine Minute lang bei etwa 200°C in der Luftatmosphäre erwärmt. Nach rascher 
Abkühlung auf Zimmertemperatur wurde die Intensität der Ph. Ex. Em. entlang der 
Proben als Abtastkurven wiederum gemessen. 


oe : 
Bei der benutzten Eintauchdauer von 5 Sekunden bewirkt Azeton nur eine verhältnissmässig 
kleine Abnahme der Erregung zur Ph, Ex. Em, 
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Auch bei diesem Verfahren wurden wieder qualitativ verschiedene Verläufe 
der Abtastkurven erhalten, je nachdem ob die Oberfläche zuerst mit einer Stahlbürste 
bearbeitet oder mit Stahl abgeschabt wurde. 

a) Mit Stahlbürste bearbeitete Oberfläche. In diesem Fall wird der Kontrast 
in den Intensitäten der Ph. Ex. Em. zwischen den eingetauchten und nicht eingetauch- 
ten Teilen der Proben auch nach der Erwärmung beibehalten. Dieses wird durch 
die Kurve I (nach dem Eintauchen in der Flüssigkeit der 2. Gruppe) und die Kurve II 
(nach darauffolgender Erwärmung) in der Abb. 4a illustriert — vergl. auch Abb. 2b. 
Das Intensitätsniveau der Ph. Ex. Em. wurde entlang der ganzen Probe durch 
die Erwärmung erniedrigt. 

b) Mit Stahl abgeschabte Oberfläche. Ein qualitativ anderes Verhalten der 
Abtastkurve der Ph. Ex. Em. wurde im Falle der Abschabung der Oberfläche mit 
Stahl (Feile oder Messerkannte) gefunden. Die Abtastkurven zeigten nun eine Umkeh- 
rung der Kontraste im Emissionsbild der Proben. Dieses ist schematisch durch die 
Kurven I und II in der Abb. 4b wiedergegeben. Die Kurve I illustriert den Verlauf 
der Ph. Ex. Em. nach dem Eintauchen in eine der Flüssigkeiten der 2. Gruppe, die 


“| 


Abb. 4. Ph. Ex. Em. Abtastkurven (schematisch) für Al-Proben, die nach einer mechanischen Bearbei- 
tung der Oberfläche in eine Flüssigkeit der 2. Gruppe eingetaucht (Kurven I) und danach an der Luft 
bis 200°C erwärmt wurde (Kurven II). 
a) Die Al-Oberfläche vorher mit Stahlbürste bearbeitet. 
b) Die Al-Oberfläche vorher mit Stahl abgeschabt. 


Kurve II illustriert den Verlauf der Abtastkurve gleichen Probe jedoch nach Aufnahme 
der Kurve I und darauffolgender Erwärmung. Das Intensitätsniveau der Ph. Ex. Em. 
wurde entlang des nichteingetauchten Teiles der Probe stärker ernied- 
rigt durch die Erwürmung als entlang des eingetauchten Teiles. 

3. Einfluss des elektrolytischen Stromes 

In weiteren Experimenten wurde der Einfluss des elektrolytischen Stromes auf 
die Triboerregte Emission orientierungshalber untersucht. Zwei Al-Proben, die vorher 
mit einer Stahlbürste oberflächlich bearbeitet wurden, wurden als Elektroden in eine 
schwache KJ-Azetonlósung eingetaucht und einer Gleichstrom-Elektrolyse unter- 
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worfen. Die Stromdichten betrugen etwa 0.1 mA/cm?. Nach einer Minute wurde die 
Elektrolyse unterbrochen, die Al-Elektroden im „chemisch reinen“ Azeton gewaschen 
und in einem Luftstrom kalt getrocknet. Danach wurden die Abtastkurven der Ph. 
Ex. Em. entlang der beiden Proben (Elektroden) gemessen. 

Es stellte sich heraus, dass an den eingetauchten Teilen der Al-Proben, die 
während der Elektrolyse als Anoden wirkten, die Erregung zur Ph. Ex. Em. immer 
stark zurückging im Verleich zur Erregung der nichteingetauchten Teile der Proben. 
Nach längerer Zeitdauer der Elektrolyse konnte sogar die Erregung zur Ph. Ex. Em. 
an der Anode völlig ausgelöscht werden, während aus dem nicht eingetauchten Teile 
der Anode immer noch beträchtliche Emission in Erscheinung trat. Dieses Auslóschen 
der Erregung an der Anode illustriert die Kurve A in Abb. 5a. Die Proben dagegen, 


A 
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Abb. 5. Ph. Ex. Em. Abtastkurven (schematisch) für Al-Proben, die während einer Elektrolyse in KJ- 
Azetonlósung als Elektroden wirkten. 

a) Al-Proben, die vorher mechanisch bearbeitet wurden und als Kathoden (Kurve K) oder als 
Anoden (Kurve 4) geschaltet wurden. 

b) Al-Proben mit natürlicher Oberflächen — Oxydschicht, die als Kathoden (Kurve K) oder als 
Anoden (Kurve A) geschaltet wurden. 


die als Kathoden wirkten, zeigten immer eine Erhóhung der Erregung zur Ph. Ex. Em. 
innerhalb der Teile, die in dem Elektrolyt eingetaucht waren, was schematisch durch 
die Kurve K in Abb. 5a wiedergegeben wurde. 

Wenn eine Al-Probe, die keine Erregung zur Ph. Ex. Em. zeigte, als Kathode 
bei einer Elektrolyse z.B. in eine KJ-Azetonelektrolytlósung eingetaucht wurde, so 
wurde diese Erregung hervorgerufen. Die Erregung beschränkte sich auf den ein- 
getauchten Bereich der Probe. Die zu einer festgelegten Erregung nötlgen Zeitdauern 
der Elektrolyse hängen von den Eigenschaften der Deckschicht auf der Al-Probe 
sowie von den Stromdichten ab. 

Die Ph. Ex. Em., die in der Azetonelektrolytlösung an der Kathode elektrolytisch 


erregt wird, zeigt ähnliche Eigenschaften wie die, die durch eine Abschabung 
mit Stahl erzielt wird, 
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Weitere Einzelheiten über elektrolytische Erregbarkeit des Aluminiums sind 
bei Lewowski (1960) zu finden. 

4. Unregelmässigkeiten des Luft-Flüssigkeit Übergangsbereiches 

Die untersuchten Al-Proben zeigten ausgeprägte Unregelmässigkeiten in den 
Abtastkurven Verläufen der Ph. Ex. Em. im Bereich der Luft-Flüssigkeit Eintauchs- 
grenze. Dieser Grenzbereich mit dem unregelmässigen Verlauf der Abtastkurve 
ist etwa 1—2 mm breit. Er wurde in dem vorliegenden Bericht übergangen, wird aber 
bei Gelegenheit nüher besprochen. 


Diskussion 


Die während der Bearbeitung der Al-Oberflüche mit einer Stahlbürste enstandene, 
sogenannte Beilby-Schicht, besteht vorwiegend aus einem Gemisch von Aluminium 
und Aluminium-Oxyd (siehe Gomer und Smith 1952). Die unterschiedliche Wirkung 
der untersuchten Flüssigkeiten auf die Erregung zur Ph. Ex. Em. von Aluminium. 
das vorher mit einer Stahlbürste bearbeitet wurde, hängt also offensichtlich mit den 
chemischen Eigenschaften der Flüssigkeiten zusammen (siehe auch Wawrzyniak 
und Sujak 1960). Beim Kontakt der untersuchten Al-Proben mit der Flüssigkeit 
kann es also zu einer Oxydation der Aluminium-Elemente und zugleich zu einer 
damit konkurierenden Reduktion der Oxyd-Elemente der frischen Beilby-Schicht 
kommen, je nach der zur Verfügung stehenden Menge der aktiven Reaktionspartner: 
Sauerstoff oder/und Wasserstoff. 

Es ist durchaus plausibel anzunehmen, dass wir dabei in erster Phase mit einer 
katalytischen Zersetzung der Flüssigkeiten an der frisch enstandenen Beilby-Schicht 
zu tun haben. Die Zersetzungsprodukte kónnen dann in zweiter Phase als aktive 
Reaktionspartner die Oxydation oder/und die Reduktion bewirken, wie das im Falle 
der Bildung von H,O, in HzO in Gegenwart von frisch gebrochenem Aluminium 
und Zink gefunden wurde (Grunberg und Wright 1952, Krause 1958, Churchill 
1939). 

Es ist bemerkungswert, dass alle untersuchten Flüssigkeiten, die den Abfall der 
Triboerregung zur Ph. Ex. Em. der Al-Oberflüche (die vorher mit einer Stahlbürste 
bearbeitet wurde) beschleunigen, in ihren Molekülen Wasserstoff- und Sauerstoffato- 
me enthalten. Die Moleküle der Flüssigkeiten aber, die den zeitlichen Abfall der 
Triboerregung zur Ph. Ex. Em. verlangsamen, enthalten keine Sauerstoffatome. 
Demnach also ist es zu erwarten, dass die Moleküle der Flüssigkeiten, die sich in 
Bezug auf atmosphärische Luft gegenüber der Triboerregung zur Ph. Ex. Em. neutral 
verhalten, weder Wasserstoff- noch Sauerstoffatome enthalten. Dies ist z.B. an Hand 
des Chlorkohlenstoffs bestätigt worden. Es ist jedoch zu denken, dass auch die Flüssig- 
keiten der 3. Gruppe (z.B. Chlorkohlenstoff) das Abklingen der Triboerregung verlang- 
samen müssen, da der Luftzutritt zu der Oberfläche der Probe während des Kontaktes 
mit der Flüssigkeit unterbrochen wird. Man darf aber die an der Oberfläche (Beilby- 
Schicht) adsorbierten Luftreste nicht vergessen. Aus diesem Grunde scheint es ver- 
ständlich zu sein, dass bei den kurzen Kontaktzeiten (5 Sekunden) kein merklicher 
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Einfluss des Chlorkohlenstoffs auf die Triboerregung zur Ph. Ex. Em. zu beobachtet 
wird. 

Die kaum merkliche Eigenschaft des Azetons, den Abfall der Triboerregung zu 
beschleunigen, scheint auf der starken Valenzbildung der Sauerstoffatome in dem 
Molekül zu beruhen (C — 0). Aus diesen Grunde, wie auch wegen der grossen Ab- 
dampfsgeschwindigkeit, kann Azeton zum Abspiilen der zur Ph. Ex. Em. erregten 
Al-Proben erfolgreich benuzt werden (Lewowski 1960). 

Wenn eine frische Oberflüche des Aluminiums geschaffen wird, so zeigen sámtli- 
che untersuchte Flüssigkeiten die Tendenz den Abfall der Triboerregung zu beschleuni- 
gen (gemessen wird eine Erniedrigung der Intensität der Ph. Ex. Em.). Dieses kommt 
vor z.B. nach einer vorherigen Abschabung der Oberflüche mit Stahl, und ebenso 
nach einer voherigen elektrolytischen oder chemischen Reduktion der Oxyddeckschicht. 
In diesen Fall wird die Tendenz der Flüssigkeiten die Triboerregung zu erniedrigen 
von uns mit der Bildung von Adsorptions- sowie Oxyddeckschichten erklärt. 
Damit findet auch die gemessene Erniedrigung der Intensität der Ph. Ex. Em. ihre 
Erklürung und es wird eine Übereinstimmung erzielt mit den früheren, gut bekannten, 
klasischen Arbeiten von Suhrmann und seiner Schule über die Rolle der Adsorptions- 
schichten in den Photoemissionserscheinungen. Zu ühnlichen Feststellungen kamen 
auch kürzlich Terenin und Willesow (1958) bei den Untersuchungen der Photoemission 
von Oxyden, sowie Patrick und Choyke (1959) bei den Untersuchungen der Elektro- 
nenemission von den p-n óbergüngen in SiC. 

Wenn also nach einer Abschabung der Al-Oberfläche der Oxydationsprocess 
durch das Eintauchen in eine Flüssigkeit der 1. Gruppe frühzeitig beschleunigt wird, 
so finden wir eine lóschende Wirkung der Flüssigkeit wieder (Abb. 3). Diese beruht, 
unserer Ansicht nach, auf der Bildung einer amorphen Oxydschicht währed des 
Kontaktes der Oberflüche mit der Flüssigkeit. 

Die Adsorptionsdeckschicht, die während des Eintauchens der Al-Probe in die 
Flüssigkeiten der 2. oder 3. Gruppe ensteht, zeigt zwar keine Oxydationswirkung 
aber sie kann die Intensität der Ph. Ex. Em. durch die Absorption von Exoelektronen 
erniedrigen. Diese Deckschicht bildet auch eine Schutzschicht gegen Luftsauerstoff 
dar. Die letzterwähnte Wirkung der Schutzschicht tritt besonders bein nachträglichen 
Oxydieren (Erwärmen) der Al-Probe sehr klar hervor. Die Oberfläche des vorher 
nicht eingetauchten Teiles der Al-Probe oxydiert nämlich schneller bei der leichten 
Erwärmung als die des eingetauchten Teiles. Dieses führt zur Umkehrung des Inten- 
sitätsbildes der Ph. Ex. Em. entlang der Al-Probe (Abb. 4a). 

Bei allen diesen Ergebnissen handelt es sich um photoelektrisch emittierte, 
negativ gelandene Teilchen und nicht um ladungslose Dämpfe, die eine Emission 
leicht vortäuschen können (Sujak 1958 a, b, 1959 b, Sujak und Bohun 1959). 

Es folgt aus den vorgelegten Ergebnissen, dass die photostimulierte Exoelektro- 
nenemission einer „frischen“ Oberfläche des Aluminiums von der zu unterscheiden ist, 
die als Folge einer Bildung einer misch-Oxydschicht (Beilby-Schicht) an der Alumi- 


niumoberfläche zum Erscheinen kommt (Sujak und Mader 1958, Böjko, Pirög und 
Sujak 1958, Sujak 1959a). " 
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Unsere Ergebnisse weisen obendrein auf die Möglichkeit der Heranziehung der 
Exoelektronenemission zur Untersuchung physikalisch-chemischer Prozesse, die sich 
an in Flüssigkeiten eingetauchten Elektroden abspielen (Lewowski 1960, Wawrzy- 
niak und Sujak 1960). 
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The contribution to the binding energy of the hydrogen molecule due to the coupling 
between nuclear and electronic motion has been calculated. The wave function employed 
was in the form of an expansion in elliptic coordinates. The most reliable result obtained 
with a four term expansion is AD’ = —5.1 cm 1. The relativistic corrections calculated 
with a five term expansion change the binding energy by AD = —2.4 — E, cm 1, where 
E,, the classical relativistic correction to the interaction between the electrons, is of the 
order of —1 cm ! or smaller. The calculated corrections added to the accurate non- 
relativistic result for infinitely heavy nuclei give the theoretical binding energy D — 
= 38280 cm 1, the experimental value being D = 38286 + 6 cm}. 


1. Introduction 


The lowest eigenvalue of the non-relativistic Hamiltonian for the H, molecule 
with fixed nuclei has been computed very accurately by Kolos and Roothaan (1960), 
but to compare their result with experiment one has to know the magnitude of the 
contributions to the binding energy due to the motion of the nuclei and to the relati- 
vistic effects. Using simple electronic wave functions the first of these corrections 
have been calculated by Van Vleck (1936) and by Dalgarno and McCarroll (1956), 
and some of the latter by Ladik (1959). The wave functions employed in these calcula- 
tions give reasonable values of the total energy of the molecule, however, the results 
obtained by Dalgarno and McCarroll show that the corrections accounting for the 
coupling between electronic and nuclear motion are very sensitive to the form of the 
wave function. It is also well known from atomic calculations that the relativistic 
corrections are more sensitive to the form of the wave function than the energy itself. 
It seems therefore desirable to calculate these effects with more accurate wave functions 


ıhan those used by previous authors. 
(129) 
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We have calculated the corrections for nuclear motion and the relativistic ones 
using for the wave function an expansion in elliptic coordinates. Only 1, 2, 3 and 4 
term expansions were employed, however, the method can easly be extended to 
longer expansions. The relativistic corrections have been computed also with a 5 term 
wave function in which the fifth term introduced some angular correlation of the 
electrons. Unfortunately, having no computer, we were not able to calculate the classical 
relativistic correction due to the interaction between the electrons. However, its magni- 
tude can be estimated, and the detailed method of numerical calculation is given in 


the Appendix. 


2. Energy Corrections Due to the Nuclear-Electronic Coupling 


It is well known that in the first approximation the coupling between nuclear 
and electronic motion changes the molecular energy by! 


E'—E—EQ— | W* H'V duds, (1) 

where V/ is the solution of the Schrödinger equation H^ = E,F with fixed nuclei, 
i n^ ; 

H= m 48 e d th A, + 2V,, V). (2) 


R denotes the relative coordinate of the nuclei, r, the position vector of the i-th electron 
with respect to the centre of R, and M is the proton mass 
The wave function ¥ we assume to be real and given by the expansion 


P => c (R) V; (5m. E2» N2 Pas Pa): (3) 


where &;, 7, are elliptic coordinates with respect to R and 9; the azimuthal angle of 
the i-th electron. 

To compute the expectation value of Ap using wave function (3) it is convenient 
to express this operator in terms of £;, 4,, p;, R, 0, (0 and © are the polar angles 
of R). After some manipulations one gets 


On oni ueniet 
dR= 3g t g BR RE [as eno ul | k 
+ A(R, 0, P, 5, M p1) + A(R, 0, D, Ea Nas Pa) + 


NEED Sn es 
" R? (£ — 3) (£i — nA E (63 — Dar + 1 (1 — xi) z] x 


SUN 9 


X [similar bracket with 1—2] + — z: | In 
(Si — 1) (L—)]*. 8g, 


a 


1) For the derivation and approximations used see Van Vleck (1936), Ta-You Wu and Bhatia (1956), 
and Dalgarno and McCarroll (1956). 
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+ cos p, 


(&-) an | 9 9 2 2 
Au ^ Vh 53. ja meg mes 
Um ee —1)-s O 8 \\ 
ED 35, * "Bom hag — 55 ) = 
(4) 
where 


E + 4? — ] ere > o 
A(R, 0, 9, È, N» 7) WI R2 pz n?) (GS 1, JE T a x n’) 2n + 


l 1 9? 2 9 9 
D. i [eom a) ża X R? (£8 — $2) le- 1) dE ln) 24x 


I 1 1, 0 coth 9 Ilə 
x|1+R — 2 — 
| E a 1 JS iss 0 Ip : sin 0 ID = A 09 ` 


a 254 kin 9  cosg O 
R*[(£& 1) (—$5]4 |^? 28 ^ sin 6 Jð 


ZE — — n2)|% si 
KG mno c > = ae > sn p eot0 2) ną Ez: 


+ cosp eot 0 Z) 2 — 


RE (E2 — 7?) 96 ' sin 0 26 


Tf we restrict ourselves to wave functions of the form 


p = > (R) 84 (415 7A > a5 No) vee (R5 7A > Es, Na p) (5) 


where p = pı — 9, and R? X? is a homogeneous function of the p-th degree of r; only, 
then the expression (4) can be reduced as follows: 
From the definition of 4? we have? 


Vr R? xy? =0, (6) 
and 
J 8 grido =>) | EREA ADR Peg, dr, dz (7) 
kp 


Now, in (7) we can drop all terms containing derivatives with respect to any angle, 
and due to the symmetry 4? (p + 2a) = 4? (p) all terms linear in cos 9, or sin g,. 
If in addition we take into account the fact that the wave function is symmetric in 
the coordinates of the two electrons, we can replace Ap by pe were 


8 209 2+- AE : 
sg + ROR? Rew ee C Vz: 


ateng- ond 
— (= 6 Dg hl M) s T+ RAR + 


Ag = ZEP maa) 


O E rM 


OE p means differentiation with respect to R with r; kept constant. 
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* RG — em EG 1 zę, t(l- m) 5 ED zę 


+ na (1 — m) z + * [(6? — 1) (1— ni) (£8 — 1) (1— )]” 


o ð 2 2 
x cos (Pı — gs) Ee mk A | ZY et 


12 dé, = Śą (8) 


Finally, to compute (7), one has to know the derivative 


> d 
= 3 Fp Pe h sms Pr Pa) 


which can be found in the following way. Let us denote 


(9) 


Bic [ PHP. Vą= | V, VT, C C VET 


(10) 
For any R the coefficients c, (R) satisfy the set of equations 
2 Ce (R) [Hix — E, 54] = 0 
k 


> Ci Cp S E IE 
i,k 


differentiation of (11) gives than after simple manipulations 


(11) 


R lS 


PNE (ENER Z 
» [Ea — B, Sa] = E -Y Hla Va. dk | 


3x do, Se = 3 
‘dR m 
i,k 

dE 2&—V 
"4,377 


SE (12) 


and by solving this set of equations one can easily find the derivative de,/dR. Simila- 
rly, by differentiation of (12), it is possible to get any higher derivative of V with 
respect to R. 


Let us now define 


, h? 
kg = — 3 PB Pr, Pr, dr dr, 


A convenient expression for the calculation of E, we shall give under the assumption 


Ba = Epp = CTT S E EUR 
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and for two special functions 


ri \2 
y? = 0? = (223) ; (13) 
and 
xy? = D = [Ei — l) (l — n (E 1) (1 — FE cos p (py — pə). (14) 


Due to the symmetry we have 
E; = P aber | "dać Vr V, En xh 2 (Wa x’) (V nEn) R En Vr, Vr, je; dt dr,, 
So, using the relations 


a 20er 
V r V v, £r = —pz (E — x?) (& — 72) 


Grsrmtm- 2 = o”) (a & — 


— k+ 0) (afa — W + 0) + 


Ag, i Sel 


+ g(B—4(B—) 


(ams £y — kni + lei) (an Ea — k' na + UES) + 
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in the case (13) we can express E, in terms of integrals 


5 "i at ns 
—2al, + 63) gs — —1 102 — dt, dT, 
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where s > —1 and k, l, k, l 72 0, and in the case (14) in terms of 
rer tni EF na dE, dE. dma dng 


The last conclusion follows from the orthogonality of Ø? functions with different p 
indices, and from the identity 


P =é t ni + 62 + 7a — 2 — 2667172 dE. 


3. The Relativistic Corrections 


For the calculation of the relativistic correction £j, we have used the Hamilto- 
nian given by Bethe and Salpeter (1957). So, in their notation 


Eq [PH QW dec Ep E, Es Ea + Ey t Ee. (15) 


For the ground state of the H, molecule E; and E, are zero, and if use is made of the 
symmetry of W the non-vanishing contributions to Æ e read 


hi ‘ 
Ei rur enn: 9 (4 n PY dt, dt; . 
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4a? e? h? 
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Eas ns f VV 6 (144) dt, dzy, (16) 


where p, is the momentum operator of the i-th eleetron, and ry, the position vector 
of the electron denoted by 1 with respect to the nucleus a. 

The calculation of E, and E; is straightforward. If the function (13) is used in 
the calculation of £y it gives integrals with 775° in the integrand. We have not been 
able to find a method of integration in this case. If we use the function (14) the integra- 


tion in Æ} can be performed in the following way. A typical integral for the computa- 
tion of FE, reads 


J Ars; 99) (A, e?) dz, dr, (17) 
As a consequence of 


A, ^ = 0, 


\ AU HEAT 3 
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after integration over g, (17) vanishes unless p = p’. Now it is clear that (17) is a sum 
of integrals of the type - 


[rang (E— 107 (EE nd 6 dt, dn, dry. 


The £,, N integration is trivial, and the integration over 74 leads to two types of 
integrals 


ee e, 
i 
&-—1 
—2a8, £s ]o 21 p 
1 
which can be expressed in terms of Ei (—2a) if s~ —1. In the case of s = —1 one 


has to perform the integration numerically. 


4. Numerical Results 
„The wave function we employed was of the form 


n 


1 A 
V gcn > cg(R)g, n=1,2,3,4,5, (18) 


where 


De EST HATTE ae, 
gą = é + £g, gy =2[(6 — 1) (1 — ma (& — 1) (1 — m)” cos (p, — 93). 


All caleulations have been carried out for the internuclear distance R = 1.4 au, 
and using « = 0.95. The coefficients c; (R), and Ej, i.e. the expectation values of Ho 
(in au), as well as the corresponding binding energies D (in eV) are given in Table I. 


Table I 
n E, (au) | D (eV) Cy ĉa CH c cz 
1 —1.08612 2.343 4.07284 
2 —1.11544 3.141 3.45343 2.03585 
3 —1.13560 3.690 3.38893 2.16629 — 1.47897 
4 — 1.15165 4.126 1.73628 2.00821 —1.38824 | 0.82506 
5 — 1.16081 4.316 1.70621 1.95257 —1.35017 | 0.84280 —0.16518 
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The corrections due to the coupling between nuclear and electronic motion 
have been computed with four wave functions (n — 1, 2, 3, 4). They are presented 
in Table H in which the first column gives the number of terms in the wave function, 
the next three give 


pan is if V An Wd, dts, 


. h2 
Es = — Be frua. Ir A,.) Pdt, dts, 
h? 
Ez = — 3M VV, V, P dt, dt, 


in au/M (= 119.47 cm-!), and the last one E’ = E, + Ey + Ey in cm’ 


Table II 
i Fi E3 E3 E^ (ems) 
1 0.9822 0.6508 0 | 195.1 
2 1.1640 0.6790 0 | 220.2 
3 1.2188 | 0.6942 0.0270 231.8 
4 0.4401 0.5831 0.0198 124.6 
Re 0 1.4519 0.0401 178.3 
Raves 0.5 0.5 0 119.5 


Although the final corrections do not show any convergence the result obtained 
with the 4 term wave function seems to be the most reliable. In this case the second 
term, ky, which is proportional to the kinetic energy of the electrons, differs only 
slightly from the accurate value 0.5875. Van Vleck (1936) employing the Wang wave 
function obtained a smaller correction E’ = 114 cm~}, and the Heitler London wave 
function with & = I in the atomic orbitals (Dalgarno and McCarroll 1956) gave a still 
smaller value. 

In the last two lines of Table II accurate values of the corrections for infinitely 
separated and coalesced nuclei are given. By comparing our 4 term result with that 
for R = oo we see that the coupling between electronic and nuclear motion gives 
a negative contribution to the binding energy AD’ = — 5.1 cm-!, whereas Van Vleck 
obtained a positive value AD’ = + 5 em-}. 

The relativistic corrections have been computed with all five wave functions (18). 


rp ^ " A d 2 
Ihe last one accounts for 92%, of the binding energy. The results in cm-1 [S au 


= 5.844 em 1) are presented in Table III. 


Table III 

n Ei E, E, | E (cmrt) 

1 — 15.3 12.2 4.2 Ti 

2 — 22.8 17.6 4.1 — 11 

3 — 24.5 18.7 3.2 — 2.6 

4 — 18.0 15:2 2.4 — 0.4 

5 — 17.7 15.0 22 — 0.5 
Ladik — 46.3 15.5 1.9 —28.9 
REQ —158.1 129.0 7.8 —21.3 
R= d — 14.6 11.7 0 — 2.9 


In the last column of the Table the sum E” = E, + E, + E; = E — E, is given. 
Accurate values for R = 0 and R = co are also shown in the Table. As mentioned 
in the Introduction we were not able to calculate Ey. This correction, being due to 
the retardation of the electromagnetic field produced by an electron, vanishes for 
R = oo and has the value E, = — 1.64 cm”! for R = 0. For R = R, it is probably 
of the order of —1 or smaller in absolute value. Ladik's (1959) results obtained with 
the Wang wave function are also given in Table III. The origin of the discrepancy 
between his result for E, and ours is not clear?. 

From Table III we see that the relativistic corrections give a negative contribution 
to the binding energy AD = —2.4 — E, cm-t. Thus, according to our calculation, 
the relativistic effects and the coupling between electronic and nuclear motion decrease 
the theoretical binding energy D, = 38287 cm-t by 6—7 cm-!. The final result is 
still in agreement with the experimental value D,,, = 38286 + 6 cm-t, however, 
more accurate calculation of the corrections is clearly desirable. 


Appendix 


Now we shall outline a method convenient for the calculation of E, with a wave 


function of the form (5) 
Pr Cu id > ch 84 (Ś1> > 52 No) PP, l (A1) 


k,p>0 


where g, are polynomials, and 
DP = y? = A? cos p (9, — Po), 


A, Um) NEE. 


3) In Ladik’s paper £j, is decomposed in a slightly different way (Ej, = E, + E,, Er = E;, 
ud z-ECRAE, 190 iL = £,,where E;, means E; in Ladik’s paper), and the relativistic Hamiltonian 
he uses is different from ours, E, being counted twice. This, however, does not influence the results 


given in Table II, 


138 


It can be shown that £j is expressible in terms of the well known W integrals (Kotani 


et al. 1955), and of the integral 


l 
i- [v1 Ud, dtTą, (A2) 
where 
2rı2 exp [— « (5 + £9)] 
Ud EIE. (Qiu, + Ogu 
R (& — 7?) (E — 7?) (2 — 73) | 1% 2 aż. 
The functions u, have the form 
u; = D) afa, E TA śą ng OP, 
and4 
O, = 511 $a tla; Os = D n en: 
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pp'kkI = 080 EP aol EF nh Dr Dr P dé, ... dy (A3) 
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DP PP — Dtr 4 4 Gre’, p>p’, (A4) 
also in terms of 
T a —2a(5,+5) ER pl ER’ yl’ 428 Oi je 
psk = | e 290092 EE yt EF ną At: OP ry d£, ... dg. (A5) 
Now we can expand o-? in a series. We have obviously 
l d 
42 SE. cS zł 
e A sin p dp ge (A6) 
where 9 = 9, — gs. 
By making use of the well known Neumann's expansion for o-! one gets 
2 lec 5 * " sin np 
03 = 7 Y nD rn Qt (£4) Pe (£-) Pr (m) Pr (a) Uwe (A7) 
t=0 n=1 
with 
Dz, = (—172(2: + 1 Es zd 
4 2 4 
) Note that O, = — (4 — zy), 0, = (ry — r$), hence (A2) is finite, 


R R? 
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The substitution of (A7) into (A5), after integration over pı an Qe, gives 
B n 
Ipskei = diae E o ICE eT] dps (A8) 
where 
DD, | O; e 18669 gh gd EE of tr Qn (E,) x 
x Pr (&_) P? (n) P7 (ng) d& ... dus, (A9) 
and j is an abbreviation for kk'll. 


Remembering the definition of O; one sees that (A9) is a combination of integrals of 
the type 


R e e MAO ET DES 
(A10) 


From (A9) it is clear that s > — 1. So if we remember that $, > £. and n >lwe 
can easily show that (A10) does not diverge 
Let us now define 


1 s 
Kr, = [ yn? Pt (n) dn. (A11) 
FE 


With this definition (A10) reads 
Er re Z (8 —1)* Ex 
x PZ ($_) de, dés, (A12) 


and making use of the generalized Ruedenberg’s formula given by Harris (1960) 
we have finally 


dE ebooks gite [ GA, (0) Gh (2) dx, (A13) 
where 
ee Ps (oy "| et 6 (po) Pr (yg 
1 


It is easy to show that 


exp (— 2a) T. 
lim Gz (x) — n-4-1 s : (A14) 
a 0 ifs — 1. 


From (A14) it follows that the product G7,, G%,y has no poles and (A13) can be simply 


computed by numerical integration. 


5) If s > n then (A10) is expressible in terms of the W > integrals. 
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The danger to aircraft of ice formations was given rise to a number of experiments 
in which the ice formation is produced artificially under controlled conditions of tempera- 
ture, airspeed, size and number of suspended water particles. 

The applications of the results of such experiments depend on information about 
the number and size distribution of water particles in natural clouds. 

In order to get such necessary information the author on the behalf of the Air Mi- 
nistry, Meteorological Research Committee, carried out during the war the investigations 
in all types of clouds over the British Isles and fogs (Fog Investigation and Dispersal 
Operation better known as FIDO). 


A 


Introduction 


Cloud is a colloidal suspension of water in air and has been described by Wigand 
and Schmauss (1929) as an aerosol. The suspended water is in the form of liquid 
droplets, ice crystals or mixture of droplets and crystals. Water particles in clouds 
are often supercooled in a liquid state at temperatures far below freezing point, some- 
times at —50° C. (Perrie, 1950). 

The danger to aircraft of ice formation has given rise to a number of experiments, 
in which the ice formation is produced artificially under controlled conditions of tem- 
perature, airspeed, size and number of suspended particles. The application of the 
results of such experiments depends on information about the number and size 
distribution of water particles in natural clouds. In 1941—42 the writer prepared for 
the Director, Meteorological Office, Air Ministry, London, a memoran dum, summarizing 
the existing knowledge on this subject and indicating the lines, on which further 
investigation was needed (Mazur, 1941, 1942). As a first step it was necessary to devise 
an apparatus for sampling the water droplets, to develop a technique and begin ob- 
servations. It is worth to mention that the first systematic study of cloud drops 
was carried out by Bricard (1939, 1940) on the Puy de Dôme and later on the Pic-du- 
-Midi, 

(141) 
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General description of the Apparatus 


The apparatus, known later as the Aircraft Impactor, was conceived on the 
lines of the Porton Impactor (Davies, 1942). The latter consisted of a wood and metal 
holder for a glass plate. The wooden front is “streamlined” and has a slot through 
which the natural wind passes, thus playing on the surface of the glass plate clamped 
at right angles to the outlet of the slot. The air stream impinging upon the plate escapes. 
round two sides of the plate. The air has to turn through a right angle, and the particles 
having sufficient kinetic energy are unable to turn the corner and so strike and adhere 
to the plate. Smaller drops will be more or less deviated but all drops greater than 
some critical size will strike the slide. 


Fig. l. Crown copyright reserved. Reproduced with the permission of the Controller, H.M.S.O 
Front view of Aircraft Impactor showing front shuter open. 


* 
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The Aircraft Impactor 


The first type of aircraft impactor, made in the Instruments Branch of the Me- 
teorological Office at Stonehouse, retained the shape of the Porton impactor but 
differed in regard to the shutters. Those consisted of two coaxial semicylindrical cases, 
arranged to rotate in opposite directions and each provided with a slot parallel to the 
axis. By suitable setting of a spring controlling the mechanism the slots could be 
brought together to give the correct exposure. 

After trials of this first model in a wind tunnel and on aircraft an improved impact- 
or was made. The body was cut from a solid aluminium block and was therefore very 
strong, the slot was shortened and the sampling plate reduced to half the size of 
a quarter plate (Mazur, 1943). 


Fig. 2. Crown copyright reserved. Reproduced with the permission of the Controller, HM.S.0. 
Back of the Aircraft Impactor showing plate holder swung round to position for changing plate. 
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Photographs of the instrument are shown in Figures 1, 2, 3. 
Figure 1 shows the front of the instrument with the front shutter open. Behind this 
shutter there is another shutter, whose time of opening can be calibrated by means 
of a cathode ray oscillograph and a tuning fork. Of the two cables shown one is a “rele- 
ase cable"; when it is pulled the shutters work in the following order: 

(1) The front shutter opens and is held by a stop. When it reaches the stop, 
the rear calibrated shutter is released. 

(2) The rear (or calibrated) shutter exposes and when the exposure is finished, 
the front shutter stop is released and the calibrated shutter is held by the recoil check 
lever. This is to prevent the shutter bouncing back and re-exposing. 


(3) The front shutter closes. 


Fig. 3. Crown copyright reserved. Reproduced with the permission of the Controller H.M.S.0 
Aircraft Impactor: view with cover removed to show the actuating mechanism 
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When the other cable — the “rewind” or reset cable is pulled, the following operations 
take place: 

l. The recoil check lever is disengaged. 

2. The calibrated shutter is rewound. 

3. The front cylindrical shutter is rewound. 

4. The rewind plate returns to its normal position when the rewind cable is 

released.. l 

The plate holder seen in Figure 2 consists of an aluminium sheet with pieces cut 
out in the vertical sides to facilitate changing plates. In the operational position this 
plate holder is mounted on ball-bearings (taken from a German plane shot down over 
Salisbury Plain) and can be turned easily when plates have to be changed. 


The Sampling Plates 


The surfacing of the plates used for the sampling of the drops is of the utmost 
importance. The plates were washed first in hot chromic-sulphuric acid solution, then 
rinsed in distilled water, allowed to dry in a dust-free atmosphere and finally polished 


R 
4 
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Fig. 4. Crown copyright reserved. Reproduced with the permission of the Controller, H.M.S.O. 


The first type of Aircraft Impactor inside aircraft. One can see the box with slides for sampling 
the drops. 
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Fig. 5. Crown copyright reserved. Reproduced with the permission of the Controller, H.M.S.O. 


The first type of Impactor on Aircraft. The window on the right side is removed. 


with untouched filter-paper. The receiving surface can then be coated with a uniform 
film of a substance in which the water particles can be trapped. Two methods of 


coating have been used: the “absolute? method and the magnesium oxide method. 


I. “Absolute” method 


In the first method the object is to preserve the drops in their original form by 
enclosing them in a surfacing substance, which must have the following pro- 
perties: 

l. a density not very different from that of water, 

2. adequate viscosity, - 

3. durability and non-solubility in water, 

4. colour contrast with water drops, 

5. non-freezing and showing very little change of viscosity with temperature. 
The matrix must be smooth and highly viscous (Benedicks and Sederholm, 1944), 


but of refractive index different from that of water. 
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The problem of proper surfacing was solved quite satisfactorily by Fuchs and 
Petrjanoff (1937) who found that good results can be obtained by collecting the drop- 
lets on glass slides coated with a freshly melted mixture of light mineral oil and vaseline. 
The proportion of the ingredients has to be varied according to air temperature and 
is 3:1 at 20° C. It is quite clear that the light mineral oil is the more important in- 
gredient; its refraction is higher than that of water. 

This method of coating was first used in Great Britain by the present writer 
during his war-time investigations on natural clouds and radiation fogs (Mazur, 
1943, see also Krauss and Smith 1949, Perrie 1950, Best 1951, Brewer 1952). The 
light mineral oil was previously saturated with distilled water and every effort made 
to obtain a bubble-free film mixture. 

The method adopted for the collecting of drops in natural clouds and radiation 
fogs (Mazur, 1944) was to place immediately after taking a sample a warmed cover- 
-glass on the underside of which was a large drop of the same substance on the layer 
over the sample; the layer melted and caused the now enveloped drops to be drawn 
up as spheres. The cover-glasses have been warmed before use in flight by putting 
them on the flat surface of a lagged vessel filled with hot water at about 45—50° C. 
When the matrix cools and sets, the droplets can be preserved for 30 hours or more 
without showing any perceptible diminution in original size or change of shape. 

In order to avoid a danger of coalescence and a piling up of droplets into layers 
when.being captured on the plates, the duration of sampling has to be very short. 
Small distance pieces were sometimes used between slide and cover-slip to avoid 
squashing the larger droplets. 

Some authors have obtained satisfactory results, using other mixtures: (i) Häuser 
and Strobl (1924) used neat's-foot oil and glycerine; (ii) Diem (1942) used two layers 
of different oils of different densities, collecting the drops on the interface between the 
layers; (iii) Bricard (1939) used mineral oil, to which organic acids had been added 
in order to reduce the surface tension to render easier penetration. Pure vaseline as 
proposed by Mason and Ludlam (1951) is definitely not suitable as a matrix for slide- 
-surfacing for the following reasons, already stated by Houghton and Radford (1938), 
see also Tahakosi (1937), Mazur (1952): the surface of the vaseline layer on the glass is 
not smooth and when viewed under dark ground illumination it is streaky and irregular, 
so that some of the captured drops are always distorted and the contrast is very poor; 
further, the density of vaseline is 0.873 g/cc at 23° C. 

The above described “absolute” method is useful as it allows the droplets to be exa- 


mined microscopically and counted at leisure. 


II. Magnesium Oxide method 


In the second adopted method MgO is “smoked” on the glass plate by means of 
a gun burning magnesium ribbon, which is moved to and fro under a glass slide at 
such a distance that the tip of the flame just clears the glass. The uneven edges of the 
layer may be trimmed by scraping off with one motion of a sharp straight edge. 
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The white magnesium powder adhers to the glass as a smooth and soft layer 
with a grain size of about 0.44. When a drop strikes the MgO layer it penetrates the 
soft surface, like a round object falling into snow, and leaves a well-defined circular 
impression; the water, of course, evaporates. Viewed by strong transmitted light the 
holes are seen as bright spots on a darker background. The MgO layer should be at least 
as thick as the diameter of the largest droplet to be measured and preferably about 
twice as thick; too thick a layer is very fragile and flakes off the glass easily, 
May 1950. 

Laboratory experiments with droplets of substance such as dibutyl phthalate, 
which can be sized accurately by other methods, enable drop size calibration; the 
ratio of diameters of drops obtained by “absolute” and magnesium oxide methods 
is about 0.85; the results are independent of drop size. 

The MgO method has advantages over the *absolute" method in that: 

(i) it is much simpler as nothing has to be done to the plate after the sample is 
obtained, 

(ii) the craters remain without change for a long period of time (2—3 days). 
Disadvantages: 

(i) the layer is very soft and fragile and can be damaged by high velocity air 
jets during exposure to the droplets; the writer found that the oxide layer is blown 
off the plate at a speed of about 220 m.p.h. 

(ii) the method is of no value for droplets smaller than 6—8u; the “absolute” 
method can show drops below 24 diameter; 

(iii) the drop of liquids other than water sometimes bounces out of the layer. 
To compare the results by the two methods the glass plates were coated half with 
MgO and half with the oil-vaseline mixture. 

The magnesium oxide method is only a modification of method developed by 
Strazhevsky (1937), who used a uniform layer of soot smoked on to a slide during his 
investigations ,on the atomization of liquid fuel. 

It is worth to mention the oldest *stain" or “absorbent paper" method, 
recently modified by Pauthenier and Brun (1940). 


Drop Counting Technique 


The droplets were counted against a Patterson-Cawood (1936) globe and circle 
comparator graticule, traverses being taken across a mechanical stage. For a good 
count about 300 drops should be taken. The number percentage in the ranges between 
the mean diameters corresponding to adjacent pairs of circles are then worked out 
and a rough curve of per cent by number up to given diameter is plotted against the dia- 
meter on a logarithmie probability paper. From a smooth curve through the plotted 
points a new number distribution is taken and converted into a table of total number 
N in the size range tabulated against size range. The mean diameter D of each size range 
is cubed and then multiplied by the number of drops in the rangeto give a table of ND3 
which, when reduced to a percentage gives a mass or a volume distribution. The 
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total mass in one hundred drops is equal to 


LC > 


M = 
6 


where o is the density in g/cc, and ZN = 100 
The concentration in g/m? can be calculated and the mean volume diameter can 


3 3 
VE Nep 
100 


The sampling errors at the very small end of the spectrum of the drop size render 


also be derived from 


impossible to take the total number of drops present in a cloud ; however the amount 
of water contained in very tiny droplets is insignificant. 

In his excellent theoretical paper Best (1951) showed that the size distributions 
of raindrops, and drop size distribution in clouds and fogs as well, may be described 


w [E 


where F = the function of liquid water in the air comprised by drops with diameter 
less than X, and a and n are constants. Best, used experimental data of the present 
writer and other workers (Diem 1942, 1946; Frith 1951, Hageman 1936). 


by the formula 


The Airspeed through the Aircraft Impactor 


In order to obtain trustworthy readings of airspeed in the aircraft impactor special 
two sets of experiments were made in the wind tunnel. For the first an old pitot 
static-tube, external diameter just under 4 inch, and spacing between static and 
pitot holes 0.65 inches, was used. Inside a slot which measures 0.6 x 1.4 inches this 
distance between pitot and static holes is obviously a serious disadvantage (Peters, 
1931). For the second set of experiments a specially fine pitot tube was constructed, 
using two stainless steel tubes each of 0.9 mm. diameter and 35 cm. length, soldered 
together. One advantage of the second type is that it does not disturb the flow of the 
air. 

Measurements with both types of pitot were made at speeds of about 70 m.p.h. 
and 50 m.p.h. The results of the calibration are shown in Figures 6 and 7 respectively. 
From the latter, in which the fine pitot tube was used, it is seen that the velocity of 
the air sample in the slot is approximately equal to the velocity of the apparatus 
through the ambient air. Many experiments carried out in wind tunnel with both 
types of impactor have shown that there is insignificant turbulence in the channel 
of the apparatus and the air velocity in the channel is very little different from the 


velocity in the undisturbed air stream. 


150 


U 


U 


Air Speea 


DD Qi AQ 
im 


S 


mites per hour 


R 
U 


= 0 7 2 3 4 5 6 7 inches 
Distance from Orifice of Jmpactor 


Fig. 6. Crown copyright reserved. Reproduced with the permission of the Controller, H.M.S.O. 
Calibration of speeds in Aircraft Impactor. 
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Fig. 7. Crown copyright reserved. Reproduced with the permission of the Controller, H.M.S.O. 
Air speeds along central axis of Aircraft Impactor (by means of a specially constructed Pitot-tube). 


Results of the observations made in Natural Clouds 


During the period January 1943 to May 1944 observations in natural clouds have 
been made over Salisbury Plain, using a Hampden bomber from Boscombe Down 
aerodrome, Royal Air Force Station. 

Consecutive samples in the cloud formations at different points both in horizontal 
and vertical positions were taken during every flight. 

The accuracy with which the drop diameter can be measured is limited by the 
resolving power of the microscope and the grain of the films. The grain size of the 
KODAK plates used is so small that it does not interfere with the definition of particles 
at least down to lu. It was possible during these experiments to see drops of about 
lu in diameter. No flattening effect of the drops was seen. 

The curves which represent the drop size distribution or the size frequency for 
a cloud are very similar to a Gaussian or normal distribution curve (Maxwell curve). 

The most obvious and characteristic feature common to nearly all of the curves 
is that there is a single and very strongly marked maximum. 

In none of the samples so far examined was there more than one peak in the size 
distribution curve. 
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The distribution curves evidently become asymptotic to the axis for larger particle 
sizes and a definite and very sharp lower limit of drop size exists. 


The clouds are not homogeneous with regard to droplet size; the diameter of 
drops vary from one part of a cloud to another. 


The locations of the maxima of the distribution curves vary not only for different 
clouds but also within the same cloud. l 

In the interior of clouds there are larger drop sizes than in the boundary layers 
where the droplets must be in equilibrium with saturated ambient air. 

The ranges of size are as follows: the diameter of the smallest drops was less than 
lu, and that of the largest above 100u. The water contents in clouds varies from 0.12 
to 0.41 g/m?. 

In well-established clouds the drops are larger and more uniform than in the 
beginning of the formation or dissipation. 


The radii of drops were found to vary in the ratio about 1:100. 


Fig. 8. Crown copyright reserved. Reproduced with the permission of the Controller, H.M.S.O. 
A micrograph of cloud drops (oil — vaseline mixture) 


Two micrographs of cloud drops embedded both in the grease film and in a magne- 
sium oxide layer are shown in Figures 8 and 9 (cumulus cloud samples). SM 
The frequency diagram for one flight, given in Figure 10, shows a typical size distribu- 


tion curve for a cloud. 
(Stratocumulus, flight on 27. 1. 1943, 2.000— 3.000 feet). 
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The numerical data of this flight are given below. 


First sample Second sample Third sample 
WARE in Number of Diameter in Number of Diameter in Number of 
microns drops microns drops microns | drops 

2 72 3 100 2 36 

256 7 188 5 120 

5 | 384 9 208 9 148 

f 300 12 176 Ti 140 

9 184 | 17 116 LĄ 88 

25 1674 24 56 26 28 

29 8 37 12 32 8 

1220 856 568 

The numerical result of the first series of flight are tabulated below. 
Diameter of the drops in microns 
Number , 
"e Genus Concentration 
of ZUeR Min. Mean. Max. 

1 Stratocumulus 2 24.2 100 0.16 g/m? 
2 Stratocumulus 2 14.4 | 37 0.12 g/m? 
3 Cumulus 9. 14.7 50 0.13 g/m? 
4 Stratocumulus 3 13 310 0.12 g/m? 
D Stratocumulus _ 15 13:3 35 0.15 g/m? 
6 Altostratus 1 10.6 60 0.13 g/m? 
7 Cumulus 1.4 8.8 21 0.12 g/m? 
8 Cumulus 3 10.2 18 0.41 g/m? 
9 Stratocumulus 4 9.0 16 0.38 g/m? 
10 Cumulus 6 12.6 18 0.20 g/m? 


The investigation was carried out at the Physies Department, Chemical Defence 
Experimental Establishment, Porton Down and all flights were made from Boscombe 
Down aerodrome. 


From the historical point of view it is worth to mention that during the 
second. world. war.investigations performed with the aid of aircraft heve been carried 
out by Professor Diem in Germany (for the Luftwaffe), by the present writer in England 
(for the Royal Air Force) and subsequently by the National Advisory Committee for 
Aeronauties and the Mount Observatory, Washington in the United States. 

(Best 1951; Diem, 1952; 1948; Krauss und Smith, 1949; Perrie, 1950). 

Some experiments concerning the behaviour of the drops, especially the problem 
of coalescence, were carried out at the Physics Department, Imperial College of 
Science and Technology, London (Swinbank and Mazur, 1943). 
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Fig. 9. Crown copyright reserved. Reproduced with the permission of the Controller, H.M.S.O. 
A micrograph of cloud drops (MgO). 
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Fig. 10. Crown copyright reserved. Reproduced with the permission of the Controller, H.M.S.O. 


Investigations of Fogs 


As well as these investigations on clouds the present writer has carried out during 
the war measurements of the number, size distribution of the droplets in different types 
of fogs by using a Cascade Impactor (British Patent No 580717) and total water content 
by using an apparatus, designed by the writer and National Physical Laboratory. 
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The need for this information has arisen in connection with the trials which have 
been carried out in the dispersing radiation fog (Fog Investigation and Dispersal 
Operation better known as FIDO). 

The investigations have been carried out during the war in the inland stations: 
Porton near Salisbury (Battery Hill, Beacon Hill, Tower Hill) — low stratus; river 
Bourne Valley, Wilts and Staines — radiation fog; coastal Royal Air Force Station 
St. Eval (Cornwall) — radiation fog, sea fog, low cloud, sea-spray. 

The results of all investigations on fogs which have been duly submitted to the 
Meteorological Research Committee, Air Ministry, London, will be published in due 


course. 
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APPENDIX 
Official opinion of the Air Ministry, London, concerning this work: 


“The results of the investigations on clouds and fogs contributed greatly to facilitate take-offs and 
landings in bad weather", Destiny Can Wait, The History of the Polish Air Force in the Second 
World War, Foreword by Marshal of the Royal Air Force, Viscount Portal of Hungerford, 
London, 1949. 

"Fog Investigation and Dispersal Operation, better known as FIDO, was prepared with the 
far-reaching collaboration of one of the Polish scientist". Destiny Can Wait, p. 341, 1949. 
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ANGULAR CORRELATION 
THEORY WITH JACOB-WICK METHOD 


By M. Micu 
Institute for Atomic Physies, Bucharest 
( Received August 5, 1960) 


The angular correlation formula for decays of the type A —> B + d, B— C + d, 
is deduced by using the development of the plane wave indicated by Jacob and Wick. 
The lack of Racah coefficients and of summations over the orbital angular momentum 
quantum numbers makes the given expression for the angular correlation simpler than 


the well-known expression (S. Devons and L. J. B. Goldforb 1957). 


1. Introduction 


Without any development in spherical harmonics some simplifications in the 
collision theory have been obtained by Jacob and Wick (1959). 

In their paper two kinds of wave functions have been used to describe the relative 
motion of two particles: plane wave functions with given helicity of two particles 
|p 4 A; >and wave functions with given eigenvalues for P, J, and for the helicity of 
the two particles |JM A445 >. 

If we denote by p the relative momentum and by 44(45) the helicity of the first 
(second) particle the matrix < JM 44, |P 444a > which relates the two wave functions 


|PAı Ag >= > [JM Ay A >< JIMA, Ap IPA, As > (1) 
JM 


apart from a factor is the matrix of the 2/ + 1-dimensional irreducible representa- 
tion of the rotation group, namely 
el 


(IM Ay hy] BA łą) = 71 Dhn- (90 — 9) (2) 


In this paper some simplifications in the angular correlation theory are obtained 
by using the same procedure as the one used by Jacob and Wick in the collision 


theory. 
(157) 
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2. Angular correlation formula 


The state of the initial nucleus A is |a; jąm4). a denote all the quantum numbers 
of this nucleus except the spin and his projection on z-axis which are denoted by j4 
and m 4. If S, (S5) is the 5 matrix which makes the first (second) transition in the cascade 


decay 
A>B+d 
B>C+d, 
the final state y is 
y = 558, la; Ją ma? (3) 


In the momentum-helicity representation for the relative motion we have 


waż 315, [ep An DE Ag |Sıl a; Ją M4) (4) 


Nw 


Since the S matrix is a Lorentz invariant the final state becomes 


Ux > (—1)'B-%8 |k, Ay Agi ky A >< ką Ap Aol S2! b;jg—Àg»(hÀ Ag |S;| a; jąMĄ) 
Ag 5 
Ec (5) 

The ket vector TR Ac; hy 4, > describes the following state 
— the free particle dy with a momentum A, and helicity A, 

— the free particle dą with a momentum ks and helicity A, 


— the nucleus C with a momentum — k, and helicity A.. 
The first momentum is referred to a frame in which the nucleus 4 is at rest while 


both the second and the third momenta are referred to a frame in which the nucleus B 
is at rest. 


The probability to find this state in the final state is 
= » C ka ką kg| Sg] bs ip — Ap X< kih Ag [Sil a: Ją may]? (6) 
B 
By using eq. (1) we obtain 


(Rh A Ag |Sı| a; jam,» — (Rh A Ag | Jamain) 4 My A | SIA la > (7) 


WA nenresenfe m PENES 5 AE : i 
514 represents, of course, the submatrix of S; belonging to a given angular momen- 
tum Jg 


If the parity is conserved in the decay A + B + d, we find 
(1, 4,|Si4l ay = MB ¢_4 4 \ Sia 
14g | S14] a) NA = Ay — Ag| 514 | a) (8) 


7] 4s Np and y, being respectively the parity factors (See M. Jacob and G. C. Wick, 1959) 
for the nuclei A, B and for the particle d}. The spin of the particle d, is denoted by Sy. 
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If the initial nucleus is unpolarized and if we do not make any observation con- 
cerning the polarization of the nucleus C we have 


pe 2 | | Digi a, (P 0 — 9)|* C A> Ac| SE | b ><A Ap 


> (9) 
For polarization insensitive detectors the eq. (9) becomes 


pA E EAM Ap-- Ac A> OF ] k } 
W = 3} (= CH Aol SEP [O> <An Aal {Al a> CRB IRA CJ aBa P, (os) 
Ack 


(10) 


For arbitrary polarization detection we introduce the efficiency matrix & and 
the density matrix o of the final state (5) and from the formula W — Spur (oe) we 
obtain 


W = BP()BLQ)OI,(p0—9) - (11) 


in which e, (£5) refers to the polarization detection of the first (second) radiation and 


Bh (1) = J; (— 17878 (4, Ag ISA a) <A, Ag ISj4l a» * C a , a lale) ôa, 


AMABAB 


—Aph—AR 
(12) 
= M (= beef CAA 1527] EE ze! „(AnleelA2> (13) 


- Ag—4 A= AC 
Aykic en 


The meaning of the angles 0, p used in the angular correlation formula is the 
following: if 6,, p, define the direction of the momentum b, in the initial (centre of 
mass) frame the angles 0, p define the direction of the momentum E, in the frame 
obtained from the initial frame through a rotation of Euler angles py, 0,, — 9, followed 
by a Lorentz transformation in the new z-axis direction. The Lorentz transformation 
is chosen in such a way that in the final frame the nucleus B is at rest. 

If the masses of the nuclei A, B and C are greater than the masses of the particles 
d, and d the Lorentz transformation is not necessary and @ represents the angle 
between the momenta of the emitted particles d, and dz. 
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PHOTO-STIMULATED EMISSION OF EXO-ELECTRONS FROM 
ANODE-OXIDIZED LAYERS OF ALUMINIUM 


By TapEUsz LEWOWSKI 
Chair of Experimental Physics, Wroclaw University 
(Received August 11, 1960; translated paper received September 25, 1960) 


It is found that photostimulated emission of exo-electrons can be obtained from 
AlO; layers by anode-oxidation of aluminium with AC current or by appropriate thermal 
processing. 

The long-wave limit of the photoemission observed lies within the visible spectral 
region. The method utilized in oxidating the Al surface was found to affect the intensity 
of electron emission. The emission is considered to be due to centres resembling the F 
centres in alkali halides. 


1. Introduction 


Kramer (1955) was the first to investigate the properties of anode-oxidized alu- 
minium layers (known as eloxal layers) by methods of exo-electron emission. He 
irradiated the A150, layers formed on the surface of aluminium with ultraviolet light 
and investigated the temperature dependence of the emission thus excited (fmding 
no such dependence). 

In an earlier note (Lewowski 1960), the present author drew attention to the fact 
that photo-stimulated exo-electron emission can be excited by processing, provided AC 
current is used in the process of forming the A150, layers. The subsequent sections 


bring further results on exo-electron emission from eloxal layers. 


2. Method of measurement 


Samples made of technical aluminium were degreased in a dilute aqueous KOH 
solution and then properly rinsed with water. Electrolytic oxidation was carried out 
in an aqueous solution of oxalic or citric acid (of about 1% concentration), utilizing DC 
or AC current. Electrolysis proceeded rapidly (in several minutes) at rather high 
current densities (of the order of 107? A/cm?), and was accompanied by electrolu- 
minescent light emission. Following oxidation, in order to remove the water, the samples 
were rinsed in acetone and dried in a current of air at room temperature (Lewowski 
and Sujak 1961); subsequently, they were placed under a point Geiger counter on 
a movable support (Fig. 1). The counter, which operated in the atmospheric air, 
was bounded from below with a metal grid (100 meshes per em?) and was provided 


(161) 
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with a water jacket circulating water from a Höppler ultrathermostat with the aim 
of keeping the counter at a constant operating temperature of 40°C (Sujak 1958). 
An accelerating voltage of 70V was applied between the sample and the grid of the 
counter. 

The method was designed so as to provide for the possibility of irradiating the 


sample with light from a tungsten bulb or a quartz mercury lamp, and of heating 


to autotransformer 


Fig. 1. Diagram of measuring circuit. 
1 — 6 V/20 W bulb; 2 — amplifier; 3 — counting rate-meter; 4— high voltage source; 5 — auto- 
matic recorder; 6 — Geiger counter with water jacket; 7 — sample; 8 — heater (kanthal sheet); 9 — 
thermocouple; 10 — copper cubes with current leads to heater and with water circulation; 11 — water 
circulation from ultrathermostat (40°C); 12 — millivoltmeter; 13 — 6 V/50 A transformer. 


the sample from 40°C to about 600°C at the approximately constant rate of about 
150°C per minute. For this latter aim, an electric heater of very small thermal inertia 
was utilized. 

The electronic apparatus consisted of high voltage source, proportional amplifier, 
and counting rate-meter with output-connected automatic recorder (Sujak 1958). 


3. DC processing of A10, 


On oxidation of the surface of the Al samples with DC current, no exo-electron 
emission at room temperature was found to occur either in the dark nor in the light 
of a tungsten bulb. The long-wave photo-effect limit for such samples lies somewhat 
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above 3000 A. If the temperature of a freshly oxidized sample, dried and placed under 
the Geiger counter is raised monotonously and the electron. emission intensity meas- 
ured simultaneously (glow-curve), the curve obtained in the dark exhibites a small 
maximum of the electron emission intensity at about 330?C. The emission intensity 
is found to increase markedly only after the temperature of 500°C has been exceeded 
(thermionic emission), as shown on curve 1 in Fig. 2a. If a second measurement is 
carried out on the same sample, the foregoing maximum tends to vanish, and electron 
emission is obtained only above 500?C (curve 2 in Fig. 2a). 


500 


400 


500 T (ec) 


Emission intensity (pulses per second ) 
a 


100 200 300 400 500 T [°c] 


Fig. 2. Glow-curves for electrons from Al,O, layers formed in the process of oxidation by DC current. 
a) | — new sample in dark; 2 — second measurement; b) 1 — new sample illuminated from bulb; 
2 — second measurement. 


The analogous temperature glow-curve obtained with simultaneous irradiation 
of the freshly-processed sample with light from the tungsten bulb exhibited two maxima 
that were situated near 330°C and 400°C (curve 1 in Fig. 2b). In a second measure- 
ment on the same sample, the maximum at 330°C vanishes, but the one at 400°C 
remains. Moreover, excitation of photostimulated exo-electron emission occurs already 
at room temperature (curve 2 in Fig. 2b). 

This excitation is especially easily accessible to observation if, originally, heating 
of the sample to 400°C and cooling to 40°C have been conducted in the dark. Thermally 


excited photostimulated emission of exo-electrons decays very rapidly (curve | in 


Fig. 3). 
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The process of thermal exeitation and the observation of electron photoemission 
can be repeated many times on one and the same sample. At a fixed temperature of 
about 400°C, the intensity of. photostimulated electron emission decays slowly, de- 


creasing by some 20% during 30 minutes. 


4. AlO formation with AC current 


As already mentioned (Lewowski 1960), on formation of Al,O, layers with AC 
current, electron emission at incandescent bulb illumination is observed to occur 
at room temperature. The intensity of the emission thus excited decreases with the time 
as shown by curve 2 in Fig. 3. The temperature glow-curve resembles curve 2 in Fig. 2b. 


5. Effect of the degree of oxidation of the Al surface 


Experiments showed that not all electrolytes utilizable for the formation of an 
oxide layer yield simultaneous excitation of the sample to exo-electron emission during 
oxidation of Al with AC current. Electrolytes exist for which no simultaneous excita- 
tion to the emission of exo-electrons is observed notwithstanding the fact that electrol- 
ysis was accompanied by electroluminescent light emission and led to the formation 
of an oxide layer. The latter group of electrolytes comprises e.g. dilute aqueous solu- 
tions of KMnO,, K,Cr,0,, CrO;, Na4B40; and of other compounds rich in oxygen. 
Moreover, addition of a small amount of one of them to solutions of citric or oxalic 
acid yields samples exhibiting electron emission of an intensity many times smaller 
than that of samples processed in a pure solution of the respective acid. 

Attempts at thermal excitation of electron emission from A150, layers obtained 


in strongly oxidizing electrolytes were also unsuccessful.! 


6. Discussion 


The exponential decay of intensity of the photoelectron emission, and the fact 
that similar half-times of decay are obtained with different excitation methods, such 
as thermal excitation (curve 3 in Fig. 3) and oxidation with AC current (the initial 
part of curve 2 of Fig. 3), seem to point to the presence of similar centres of emission 
in either case. 

The process of electrolytic oxidation of aluminium, when conducted rapidly in 
an acid medium at rather high current densities, provided conditions for incomplete 
oxidation of aluminium and hence for the arising of Schottky type defects (Grunberg 
and Wright 1955, and Sujak 1959). The empty oxygen lattice points ean subsequently 
be occupied by electrons e.g. if the sample is subjected to electric discharges in a gas, 


1 The addition of even a small amount of Na,B,O, to the oxalic acid changes the spectral 
composition of the electroluminescent light, released during the formation of Al,O,-layer. The 
intensity of the yellow part decreases namely when compared with the intensity of the blue part of 
the light emitted. This shows that there exists a correlation between the ability of excitation of 
exoelectron emission and the spectral. composition of electroluminescence light in the case of ALOJ 


165 


or during the flow of electrolysis current in the direetion of transmission (in the 
process of oxidation with AC current) (Lewowski 1960), or in raising the temperature 
of the sample. In the last case, the electrons most probably pass from the conduction 
band of the metal through the potential barrier at the metal-oxide contact to the trap 
levels of the oxide. (The conduction band of the oxide lies lower and, if there are no 
other full levels in the oxide, the only way of explaining the process of thermal excita- 
tion of the oxide to photo-emission within the visible range is by invoking the tunnel 
effect). The trap levels thus invested, whose structure resembles that of the F 
centres in alkali halides, can become on illumination of the sample the source of 
photoelectrons presenting an effective work function lesser than that for the exit 
of electrons from Al or from Al,O, containing no such levels. 

Electron emission of longer time of decay (curve 2 of Fig. 3) seems to be related 
to other emission centres in Al;O,, or to the process of etching Al (the chemical reac- 
tion, cf. Wawrzyniak and Sujak, 1960) by citric or oxalic acid. This latter hypothesis 


7000 
800 


600 
IF 52 RING 
400 


12 24 36 48 60 t [sec] 


Fig. 3. Curves of decay of electron emission excited from Al,O, layers. 1 — thermal excitation, 2 — exci- 
tation by formation with AC current. 


seems to be the more probable one, as photo-electron emission of very great intensity 
and a similar time of decay was found to occur e.g. on immersing an Al sample in 
a dilute solution of orthophosphoric acid, H;PO,. On the other hand, electron emission 
occurring also in the dark at about 330°C in new samples is presumably related to 
the presence of water in the oxide lattice, notwithstanding the fact that it takes place 
only in the presence of an accelarating potential between the sample and the grid 
of the counter. Further investigation is proceeding. 

The author thanks Professor B. Sujak for his discussion of the subject and for 


his valuable advice. 
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THERMAL CONDUCTIVITY CHANGES 
OF ZnS0, : 7H,0 BELOW 10°K 


BY B SUME 
The Clarendon Laboratory, Oxford, England 
(Received Semptember 5, 1960) 


Preliminary results have shown, that treatments with water and the presence of 
acetone vapor in the gas atmosphere in which the single crystal specimens of ZnSO,-7H,O 
were stored at room temperature, influenced their thermal conductivity measured at low 
temperatures below 10°K. The water increased but acetone vapor reduced the thermal 
conductivity below that temperature. The thermal conductivity values above 10°K always 
remained unafected. 

Reductions of the initial thermal conductivity of about 5095 due to acetone vapor 
could be observed at 3.5°K. 

The observed phenomenon is thought to be due to changes of perfection of the 
specimens used, as a result of adsorption and diffusion of acetone and water molecules 
from and into the specimens which lead to changes in stechiometric composition of the 


specimens. 


Introduction 


Bijl (1949) found in a course of experiments on Chromium Potassium Alum 
[KCr(SO4),12H,0] that the curves of the thermal resistance versus temperature 
resulting from different runs, but obtained with one rod did not coincide in the temper- 
ature region below 4^K. Higher values for the thermal resistance were found when 
the crystal was cooled down slowly. Bijl (1949) concluded that the observed phenom- 
ena were characteristic of the substance examined, rather than a consequence of 
the imperfect physical purity of the rods used. He further stated that it seemed to 
be feasible that the water molecules in the crystal lattice were responsible for this 
behaviour. 

In this report some qualitative, preliminary results concerning similar behaviour 
in Zine Sulphate Hydrate will be discussed. Because the behaviour of Zinc Sulphate 
Hydrate seems to be due to the changes in the physical purity of the samples used 
(especially the water molecule distribution), some general theoretical consideration 
on the thermal conductivity problem will be first outlined. 


* On temporary leave from Zaklad Niskich Temperatur Instytutu Fizyki Polskiej Akademii Nauk, 


Wroclaw. 
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Theoretical considerations 


In the total lattice thermal conductivity K, the effects of transport of heat energy 
by waves of all different lattice frequencies are observed. To obtain a simple expression 
for the thermal conductivity it is necessary to average the frequency over the whole 
spectrum. It is convenient to introduce the kinetic formula 


R obl (1) 


Where c is the specific heat at constant volume, v is the average phonon velocity and A 
is the phonon mean free path. By the mean free path of a phonon, we mean the 
distance a wave travels before its intensity is reduced to l/e of its initial values (Zi- 
man 1960). 

At the lowest temperatures, below the temperature of the thermal conductivity 
maximum the mean free path of the phonons is mainly determined by scattering of 
phonons by grain boundaries and the external boundaries of the specimen and lattice 
defects as impurity atoms or vacancies and dislocations. 


At higher temperatures the mean free path of the phonons is almost entirely 
determined by the phonon-phonon scattering (Peierls 1929, 1935) and the recently 
measured spin-phonon scattering (Rosenberg and Sujak 1960, Orbach 1960). 

When the possibility of interchanges with the water of crystallization structure 
in a hydrate is taken into consideration it seems reasonable to assume that the changes 
in the specimen will take place first in its surface layer. This means that point defects 
(vacances or impurities) and different grain boundaries may be created in this layer. 
By extracting some of the water molecules from the specimen (hydrate) some vacancies 
will be created. By diffusion of different vapors different impurity defects will be created 
in the studied specimens. By these treatments the surface roughness will also be changed 
and this will alter the surface reflectivity. 

It may be shown (Ziman 1960) that, for phonons, the surface reflectivity may 


influence their mean free path A, as follows 


re 
een (2) 
1—p 

where p may be defined as the polish of the surface and gives a measure of the reflectiv- 
ity of the surface. For a perfectly rough surface p = 0 and for a perfectly smooth 
surface p = 1 is put. 

On the other hand if A; is the mean free path associated mainly with the internal 
scattering by lattice imperfections, Ziman shows that in addition to boundary scattering 


(by assuming that the walls of the specimen are perfectly rough) the total mean free 
path of the phonons is given by: 


3 — 
A= gis ff Atmen (- tan as Je Kin 
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where 5, is the cross section of the specimen (rod), dÓ is an element of solid angle 
giving the phonon path direction (r—rg) which is supposed to make an angle 0 with 
the thermal gradient along the axis of the sample (rod). 

The following formula can be also put down as a case of Matthiesen's rule: 


Lgresobh vieil 
dint AA A, (9 


and assuming this to be correct! for the first crude approach A will be given by 


A; A, 


Putting the total mean free path A into (1) we may write the kinetic formula 
for thermal conductivity as follows: | 


migowy ey; (6) 


Any change in the lattice perfection even if it is only in the surface layer will 
influence the thermal conductivity K, because A, will be affected. One has to bear 
in mind that the phonon reflectivity of the surface layer is not the same as the reflectiv- 
ity of light at the specimen surface. Therefore adsorption, absorption and diffusion 
processes of various vapors on the surface and into the specimen during the time of 
storage at room temperature before the cooling down process starts, will influence 
the thermal conductivity values in the very sensitive temperature region below the lattice 
thermal conductivity maximum. This will occur in a way dependent on the occurrence 
of these processes for the given samples. 

The adsorption and diffusion processes in hydrates at room temperatures occur 
easily. In particular the diffusion of the water molecules occurs in both directions, 
into or out of the specimen dependent on the atmosphere surrounding it. Therefore 
it seems that the hydrates are particulary sensitive specimens for this type of investi- 


gation. 


Experimental 


a) The specimens 

The experiments were carried on with ZnSO, : 7H,0 single crystals with FeSO, - 
- 7 H,O impurity added to the solution in which the crystals were grown. The specimen 
I included about 0.25%, weight and specimen II about 1%, weight of this impurity. 
The basic materials used was of the *Analar" purity grade (maximum impurity 
0.006%,). The specimens were cut from larger single crystals along the edges to the 
size of about 3 cm in lenght and 0.4 cm in square. The thermometer contacts as well 


1 Ziman (1960) discusses the reasons for which this formula must give spurious results by a deeper 


approach, 
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as the “cold contact” were made of copper wire and were cemented to the samples 
with “Araldite”. On the “hot end” of the speeimens a non-inductive heater coil of 


about 100€) was cemented. 


b) The method and apparatus 

The method of thermal conductivity measurement as well as the apparatus used 
was the very same as already described by Rosenberg (1955), they will be therefore 
only briefly outlined now. 

Thermal conductivity measurements were taken in the temperature range between 
2.5°K and 23°K using the stationary method. One end of the specimen was attached 
to the high pressure chamber of a Simon expansion helium liquefier with the aid of 
a copper wire connection. This provided the *cold contact". The other end of the 
specimen was heated electrically. The temperature gradient along the specimen was 
measured by two helium gas thermometers using a differential manometer. The two 
helium gas thermometers were attached to the copper contacts on the specimen. 
In view of the relativly high thermal conductivity of the specimens the power input 
to the heater was of the order of 10 to 20 milliwatts, which gave a temperature difference 
which in the most cases was below 0.1? K. In a few case it was as high as 0.4? K. 

The measurements were carried out on two specimens. The cooling down process 
from room temperature proceeded in the following way. The atmospheric air surround- 
ing the specimen mounted in the apparatus at room temperature was first pumped 
after the apparatus (and the sample) have been cooled down to about 90° K. By this 
one should avoid dehydrating the specimen. The air acted also as a heat exchange 
medium (a very poor one) in the cooling down process to about 90? K. The time 
needed to reach 90? K from room temperature was about 4 hours. After that the air 
was evacuated and replaced with helium exchange gas at a small pressure to continue 
the slow cooling down process. Liquid hydrogen temperature (20? K) was reached 
by the high pressure chamber of the expansion liquefier in about 1 hour. At the 
temperature of liquid hydrogen the measurements began. They lasted always about 
4 hours until the helium bath temperature of about 2? K was reached. 

c) The results 

The results obtained with the specimen I (0.25%, Fe impurity) for different, 

subsequend runs are collected in Fig. 1. 


Rum I 
The curve denoted by I (+) gives the thermal conductivity results for the first run 


Run II 


After the measurements of the first run were ended the apparatus was warmed 
up to about 90° K. This took about one hour. The atmospheric air was let into the 
specimen space to a pressure slightly below one atmosphere. Then the temperature 
of the apparatus was raised to room temperature and the pressure of air in the specimen 
chamber was continuously aqualised to atmospheric pressure. After this the apparatus 
was cooled down again to liquid oxygen temperature (90° K). 
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In about 24 hours the air was pumped from the specimen chamber and the helium 
exchange gas admitted. The apparatus was cooled down to liquid hydrogen temperature 


(in about 1 hour) and the measurements of the thermal conductivity started for the 
second time. 
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Fig. 1. Measured thermal conductivity k versus temperature T for sample I — ZnSO4.7H,O (0.25% 
FeSO,) (+) run I, (-) run II, (x) run III, (Bi) run IV, (o) run V (for details see text). 


The results obtained are shown as curve II. It can be seen that the values of the 
thermal conductivity are higher in the region below about 9° K than in the case of 
the first run. The maximal increase of about 10% is to be seen at about 5° K. 


Run III 

After the measurements of the thermal conductivity values given by curve II 
were finished the apparatus was allowed to warm up to room temperature over-night. 
Then the specimen was wetted with distilled water and quickly dried with absorbent 
paper. After this treatment the surface of the specimen became very clear. The spec- 
imen chamber was resealed and the cooling down process started. The results obtained 
in this third run are given by the curve III (x). 

The thermal conductivity showed a further increase especially at the lowest of 
measured temperatures. The increase is of about 30% at 3°K (the values given by 
curve II are taken as 100%). 

Run IV 

After the measurements of the third run were taken the apparatus was allowed to 
warm up again over-night to the room temperature and the specimen chamber reopen- 
ed. A small amount of acetone was put in the specimen chamber (it did not touch the 
specimen as a liquid) and the apparatus was prepared for the next run, 
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After about 5 hours the cooling down process could start. The results of the 
following thermal conductivity measurements taken are given by curve IV (m). 
A very marked decrease of the thermal conductivity below 10°K appeared with 


respect to the values of the curve given by run IH. 


Run V 

The apparatus was warmed up again to about 200°K and dry air at atmospheric 
pressure was then admitted into the specimen chamber. After that, the apparatus 
was warmed up to room temperature and the pressure in the specimen chamber 
equalized to the atmospheric pressure again. Afterwards the apparatus was cooled down 
tp 90°K (this took about 4 hours). It was kept at 90°K for 48 hours, when the next 
run took place. 

The results of measured thermal conductivity in that run are given by the curve 
V (o). 

As it is to be seen a further decrease in the conductivity accured. The decrease 
in the conductivity was more pronounced at 3^K than at 5^K just the opposite of 
the case of run IV (curve IV). 

It is very significant that all runs gave very coincident results 
from 20°K to temperatures down to about 10°K. 


High vacuum conditions 


To prove that in high vacuum conditions and at sufficient low temperatures 
(say 90°K) the specimen can be stored and the subsequent results will coincide 
completely the next experiment was performed and the results are given in Fig. 2. 
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Fig. 2. Measured thermal conduetivity k versus temperature T for sample II — ZnSO, 47H30 1% FeSO,). 
() run I, (x) run I (for details see text). 
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The specimen II (1%, Fe impurity) was used for this purpose. The thermal 
conductivity was measured in the same conditions as run I in the measurements 
described above — curve I in Fig. 2. Afterwards the apparatus was warmed up to 
90°K and keept at that temperature for about 30 hours. During that time the high 
vacuum was maintained in the specimen chamber (order of 10-5 mm Hg). 

Then the second run took place. The thermal conductivity curve obtained 
(curve II) coincided now with the former curve I even below 3°K. 

Further runs proved also that always when the specimen lost some water of 
crystallization (for example after treatment with warm air) the surface became less 
clear and the conductivity curve was reduced on the low temperature side of the 
maximum. 


Discussion 


The qualitative results show clearly that in future studies of thermal conductivity 
of dielectrics (and especially hydrates) at very low temperatures below the conductivity 
maximum the history of the specimen before the measurements take place is essential. 
This concerns not only the crystallization process but. even the conditions of the 
atmosphere at which the samples are stored. 

The changes observed in the thermal conductivity values for the same sample 
may be ascribed to the changes in the water molecule-concentration in the samples 
which may lead to changes in crystalline structure or the similar behaviour, i.e. the 
creation of point defects by irradiation or the diffusion of alkali metals into the alkali 
halides — the so called F centers (Berman et al. 1955, 1959; Pohl 1960). 

Further study is essential to find out how fast the sensitivity of the thermal 
conductivity measurements in the region below the maximum can be evaluated for 
adsorption, absorption and diffusion problems in dielectric crystals, especially hydrated 
crystals. 

The author wish to thank very cordially Dr H.M. Rosenberg for the patient and 
essential guidance in introduction into the basis of the technique of thermal and 
electrical conductivity measurements at low temperatures and for the provision of 
the samples. 

Opportunity is also taken to express gratitude to Professor B. Bleaney for his 
kind wellcome and permission to work at The Clarendon Laboratory. 

The stay at The Clarendon Laboratory was financed by the Polish Academy of 
Sciences for which best thanks should be expressed. 
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The probability of non-excitation of atomic electrons in the beta-decay of He$ has 
been calculated using approximate wave functions including the self-consistent field 
ones. The results are significantly different from those obtained with wave functions which 


give more accurate values of the energy. 


The atomic excitation and ionization which accompany f-decay of the nucleus 
have been studied theoretically by many authors. However, the mean f-imduced 
ionizations estimated theoretically do not seem to agree with the experimental values 
(see e.g. Schwartz 1953). In most of these calculations the atomic wave functions 
employed were constructed from one-electron hydrogen-like orbitals and it has been 
suggested that the more accurate self-consistent field (SCF) wave functions should 
be used to improve the theoretical results. 

The purpose of the present note is to show how sensitive are the probabilities of 
the f-induced transitions to the form of the wave function employed in the calculation. 
The most convenient quantity that can be used to test various approximate wave 
functions is the probability of finding in the ground state the Lit ion resulting from 
the ß-decay of He9. For both systems many approximate wave functions can be found 
in literature. 

It has been shown by Feinberg (1941) that in ß-decay the main effect comes from 
the shaking of the atomic core due to the change of nuclear charge. Thus, in the 
sudden perturbation theory, the probability of finding the daughter Li* ion in its n-th 
quantum state (in the discrete or continuous spectrum) is given by 


Pon = | Ii Po V, drl’, (1) 


where gy is the ground state wave function of the helium atom, and y, is the wave 
function for the n-th state of the Lit ion. Py, obviously vanishes if pọ and y, have 


(175) 
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different symmetries. Py is the probability of non-excitation, i.e. the probability of 
finding the Lit ion in its ground state. The quantities Po, fulfill the following condition 


SI pecu (2) 


where I) denotes summation over all states of the discrete spectrum and integration 
n 
over continuum. 


For the fi-decay of He® the probabilities Pg, have been computed by Winther 
(1952). To calculate Py, he employed for He the 6 term and for Li+ the 7 term Hylleraas’ 
wave functions, i.e. expansions in terms of ry + rg, Ty — Tą and Ty = Ir, — A We 
want, however, to test simpler types of wave functions, especially those obtained by 
the SCF method, which are believed to be accurate and which are available or can 
be obtained for more complex systems. Therefore we have calculated Pog using for 
He the following wave functions 


Yor = Ny e ntn (3) 
Poz = Na x(r) x(r3). (4) 
Pos = N, [ettr JE eben tira) (5) 
Por = Net (1 + crys), (6) 
qos = Ng e 0172 [1 + e ra + eg (ry — rg)?], (7) 


where 4 represents the SCF orbitals computed by Weiss (1958) in the form of expans- 
ions in terms of the Laguerre polynomials 
1 (0) pet Oar, rasen (8) 
m=1 

For the Lit ion the functions yg. ... yg, hawe been employed. They have the 
form of gg, ... gg, respectively, of course, with different parameters which in all 
cases have been obtained by minimizing the energy. In the SCF orbitals (8) 5 term 
expansions have been used; for helium also a 3 term expansion, and in this case the 
corresponding wave function will be denoted by gg. The 6 and 7 term wave functions 
used by Winther (1952) for He and Lit respectively we shall denote by gs, and 
Yow * 

The calculation of P is straightforward and the results are presented in Table I 
where in the last column the value obtained by Winther (1952) is also given. In the 
last line we give the sums of the energies of both systems obtained with the correspond- 
ing wave functions. 

Table I 


Values of Poo calculated with various wave functions 


Wave functions | Por» Yor | Poa» Yor | Poz» Poz | Pos» Vos | Pos» Poa | Pos» Yoa | Pow Pow 


Poo 0.7247 0.7317 0.7314 | 0.7110 0.7174 0.6859 0.6698 


— Energy (au) 10.0700 | 10.0980 | 10.0981 | 10.1247 10.1516 į 10.1707 10.1827 
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The results shown in Table I are rather surprising. They suggest that the SCF 
wave functions give the worst values of the probability P. Even the wave functions 
(3), which in fact are also SCF wave functions only poor ones, give for Py a value 
that seems to be more accurate. This is in contrast with calculations of other quantities. 
Thus, it is well known that the SCF wave functions yield good results when applied 
to caleulate expectation values of the operators which depend only on the overall 
electron density distribution and not explicitly on rys. Therefore, to test in a different 
way the wave functions employed in this calculation, we computed for helium the 
expectation values of r and r? using the wave functions gy, ... Pos- The results are 


shown in Table II. 


Table II 
Expectation values of r and r? for the helium atom 
Wave funetion Po1 Por | Poe Pos | Poa Pos 
| 
r (au) 0.8889 0.9265 | 0.9273 0.9379 | 0.8970 0.9237 
T 0.9563 0.9968 0.9977 1.0090 | 0.9650 0.9937 
r? (au) 1.0535 1.1803 1.1849 1.2473 1.0772 1.1736 
tlre. 0.8827 0.9889 0.9928 1.0450 | 0.9025 0.9833 


The ratios of the computed expectation values to the accurate ones (Kolos, Roothaan 
and Sack 1960) are also given. In this case the SCF wave functions yield very good 
results, more accurate than those obtained with other approximate wave functions 
which are better for calculation of the energy. 

We must admit that probably the value of Po) obtained with the SCF wave funct- 
ions is the least accurate one. However this is not certain. A wave function with most 
variational parameters which partly accounts for electron correlation and gives the 
best energy does not always describe in the best way all properties of the system. 
Let us consider e.g. the probability of finding both electrons of the helium atom close 
to the nucleus. This probability is given by |gg (0, 0)|? dr where gy (0, 0) is the number 
obtained from the wave function gy (rj, 72) by putting ry = rg = 0. The value of the 
wave function for r = 0 plays an important role in the theory of K- or L-capture 
(see e.g. Odiot 1958). The values of gy, (0, 0) and wp; (0, 0) obtained from some wave 
functions used in the previous calculation are given in Table III, where 7^; denotes 
Po; OF Yo; for He and Lit respectively. For helium they show the same type of regularity 
as the previously computed probabilities Po). The results obtained with the simplest 
wave functions of type (3) are intermediate between the values corresponding to the 
SCF wave functions and those used by Winther (1952). 


Table III 
Values of the He and Lit ground state wave functions for rı = rg = 0 
Wave function | Por | Yoo | Wo | Ps | Ua 
He 1.530 1.783 1.800 1.418 1.388 
Lit 6.179 6.843 5.926 6.219 
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Table IV 


Values of the Kinoshita wave functions for r; = ry = 0 


UNumberof time |^ 107] PO SO YN "IS c E 
Po (0,0) | 1530 | 1388 | 13m | 1923 | 19283 | 1.9303 | 1.9309 


It seems natural to expect that when the electron correlation is taken into account 
in the wave function the probability that both electrons are simultaneously at the 
nucleus decreases. However, we see that this is not true for Li*. In Table IV we give some 
additional values of pọ (0, 0) obtained from the more accurate wave functions computed 
by Kinoshita (1957) for the ground state of the helium atom. The 6 term wave function 
used by Winther also belongs to this series, as well as the simplest wave function (3), 
and for comparison the corresponding results are included in the Table. The results 
are again rather unexpected ones. We see that more than 20 terms are needed to 
obtain a value which is probably (but not necessarily) the most reliable, and that this 
value is higher than that obtained with the SCF wave functions. Thus the effect of elec- 
tron correlation in the He atom on the ratio of the probabilities of K- and L-captures, 
found by Odiot (1958), may be reduced if more accurate wave functions are employed 
in the calculation. The computation of the probability Po should also be carried out 
with much more accurate wave functions before concluding about the applicability 
of the SCF wave functions to this problem. 

My thanks are due to Mr. A. Weiss for sending the unpublished results of his 
self-consistent field. calculations. 
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REVIEWS OF BOOKS 


BEITRÄGE ZUR PHYSIK UND CHEMIE DES 20 JAHRHUNDERTS 


Lise MEITNER, Orro HAHN, Max von LAUE 
zum 80 Geburtstag 


Herausgegeben von O. r. Frisch. F. A. Paneth, F. Laves, P. Rosbaud. Friedrich Vieweg& Sohn, Braun- 
schweig, 1959, 285 S. 


Dieses Buch ist ein Gedenkbuch, dem 80 Gebrurtstag drei Nobelpreistriger, 
Lise Meitner, Otto Hahn und Max von Laue gewidmet. Es enthilt drei Serien von 
Artikeln, deren Inhalt mit den wissenschaftlichen Forschungen der oben genannten 
Gelehrten eng verbunden ist. 

Manche Artikel berichten iiber die Geschichte der Forschung der vielfiltigen 
Gebiete, in denen die Jubilare wissenschaftlich gearbeitet haben: oft ist die Geschichte 
bis zu den neuesten Entdecungen geführt. Interessant sind auch die Artikel, welche 
die persónlichen Erinnerungen der Mitarbeiter der Jubilare beschreiben. 

Der erste Teil des Buches, Lise Meitner gewidmet, enthält 7 Artikel über manche 
Probleme der « und f-Radioaktivitát. Im ersten Artikel dieser Serie beschreibt 
A. Przibram seine persönliche Erinnerungen an das alte Physikalische Institut in 
Wien. Dann folgt ein Artikel von S. Rosenblum über die magnetischen Spektren 
der Teilchen. W. Heitler berichtet über die Geschichte des Durchgangs der y-Strah- 
lung durch die Materie vom theoretischen Standpunkte aus und endet seinen Artikel 
mit den Bemerkungen über nicht-lokale Theorie. 

W. Gentner hat die Geschichte der frühen Forschungen über y-Strahlung dar- 
gestellt. Dann folgt ein längerer und besonders interesantes Artikel von C. S. Wu 
über die Geschichte der Forschung des ß-Zerfalls, der mit der Darstellung der Nicht- 
erhaltung der Parität in schwachen Wechselwirkungen endet. 

H. Olsen und H. Wergeland berichten über Probleme, die mit Bremsstrahlung 
verbunden sind. Der erste Teil endet mit dem Artikel vom A. Flamenfeld über die 
Geschiche der Isomerie der Atomkerne. 

Der zweite Teil des Buches besteht aus 9 Artikeln. über chemische Probleme 
der Radioaktivität, die mit Arbeiten von Otto Hahn verbunden sind. Der erste Artikel 
dieser Serie enthält die Erinnerungen K. E. Zimen über Kaiser Wilhelm Institut 


für Chemie. 
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Dann folgen: ein Artikel von A. C. Pappas über radiochemische Eigenschaften 
der kurzlebigen Produkte der Spaltung der Atomkerne und kurze Bemerkungen 
von B. Karlik über Szintillationsspektrometrie. 

G. T. Seaborg beschreibt die Geschichte der frühen radiochemische Forschungen 
von Plutonium und G. von Hevesy berichtet über die Anwendung des radioaktiven 
Eisens in der Biologie. Es folgen die Artikel über die Anwendung der radioaktiven 
Isotope bei der Forschung des Kristallwachstums (N. Riehl), in der Geochemie - 
(H. J. Born) und Technologie der Farbstoffe (H. Baddenhausen, H. Gótte und 
L. Wiesner). 

Der dritte Teil des Buches ist Max von Laue gewidmet und enthält 10 Artikel 
über die Anwendung der Róntgenstrahlen in der Kristalographie und Chemie. Nach 
den kurzen Artikeln von P. P. Ewald und W. L. Bragg, die den historischen Charakter 
haben, folgt ein längerer Artikel von S. Laves über Kristallstruktur und Kristallchemie 
von Elementen und metallischen Verbindungen. 

In zwei Artikeln besprechen A. Guinier und H. Jagodzinski die Struktur der 
Kristalle und Mischkristalle. K. Lonsdale schreibt über die schwingenden Atome in 
Kristallen. Zwei Artikel von W. Hoppe und R. Hosenmann besprechen die Probleme 
der Forschung der Struktur von Kristalle und Molekülen mit Hilfe von Róntgen- 
strahlen. Diese Serie schliessen zwei Artikel von J. M. Bijvoet und G. Bormann über 
den Durchgang der Röntgenstrahlen. 

Die Mannigfaltigkeit der Themen die in diesem Gedenkbuch dargestellt sind, 
ist so verschieden, dass sowohl ein Physiker, wie auch ein Chemiker und Kristalograph 
wird sicher einen Artikel finden, der ihn interesieren wird. 


Bronistaw Sredniawa 
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MODIFIED DERIVATION OF DISPERSION LAW FOR LONG SPIN WAVES 
IN ANTIFERROMAGNETIC REGULAR SUPERSTRUCTURES 


Bv HENRYK COFTA 


Institute of Theoretical Physics, A. Mickiewicz University, Poznań 
and 


Ferromagnetics Laboratory, Institute of Physies, Polish Academy of Science, Poznań 


(Received November 21, 1960) 


In one of our previous papers (Cofta 1959) we have derived by the semiclassical 
method a general dispersion formula for spin waves in ferri- and antiferrimagnets, 
which holds for all translational (i.e. Bravais) magnetic lattices with regular! super- 
structures and which goes over into formula for antiferromagnets in the case of equal 
opposite spins. lts validity for arbitrary propagation vector k is rather apparent; 
the spin wave theory in the harmonic oscillator approximation is valid only for long 
spin waves for which the interactions between spin waves may be neglected. This 
means, that one must always confine himself to small values of k and consequently 
one use practically the approximate power series expansion for small k in all applica- 
tions of dispersion law. 

The purpose of the present paper is to show, that an direct derivation of such 
approximate dispersion formula is much simpler than the derivation of the apparently 
strict one and then its expansion. The direct way used here consists in the well known 
treating of atomic spins S,, as distributed in a continuous manner over the crystal?. 
Such an approach permits to express aproximately the spin S(m + r) in the position 
m --T by the first three terms of a Taylor expansion as follows 
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1 The rigorous definition of the regular superstructures may be found in the paper by Cofta (1960 b), 
2 See for example Nordheim (1934), page 857 (the case of ferromagnetic chain) or H. Kaplan 


(1952) (the case of ferrimagnetic spinel) 
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where the symbol & is used for dyadic multiplication of two vectors, the symbol © 
however denotes a multiplication of two tensors, which contracts the product to 
a vector (note that d?S/dm? is a tensor of the third order). 

Let us denote by A and B the dimensionless vectors of atomic spin in the one 
and in the other sublattice respectively. The Heisenberg-Dirac exchange Hamilto- 
nian leads to the following equations of motion (Herring, Kittel 1951) in the semiclassi- 


cal approximation?: 
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Here H, is the effective field of magnetic anisotropy, y — the magnetomechanical 
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ratio (y<0), m and r — the positions of lattice points, and J, the exchange (or 
superexchange) integral between spins localized in points m and m + r (which does 
not depend upon m). Hence, in the antiferromagnetic case we deal with two inter- 
penetrating spin fields A (m) and B(m). Using the expansion (1) for both A and 
B 
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mcr 
we obtain for (2): 


These equations are satished by the wave-precessional solutions: 
A = R,(j, — ij,) exp (km + ot) + Sja) 
B = R,(j, — ij,) exp (km + wt) — Sj, | S 


assuming that the field HZ, has the z-direction. The substitution of (4) into (3) 
yields two simultaneous equations for the amplitudes R, and R,: 


R,[S(F— x) + E+ ho] + R,S(G — x) = d 
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where following abbreytatione Me Used! 
E,—hyH* =hly| H, | 
F(k) -kok: 5 J,(r © r) = > J (kr)? 
G(k) =k o k: N Jr o v) => J (kr)? 
x = 2X], | 


3 As to the details see the paper by Cofta 1959 and 1960 a 
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The secular equation of (5) leads to the dispersion relation for spin waves: 


|| = VIE) + S(F — x)]? — S*(G — x)? (7) 


which holds for arbitrary regular superstructure and can be obtained also by expanding 
all functions cos (Kr) in the Ziman’s (1952) dispersion law (derived using the Holstein- 
Primakoff formalism) or in that by Cofta (1959) (derived using the semiclassical 
formalism). 

We may bring the derived formula into a more convenient form taking into 
account that for small k 


F~C<x (8) 
Then, we may neglect the k* terms obtaining finally 
|ho| = [aF(k) + bG(k) + c] (9) 
where 
a = 2S(E, — Sx) 
b= 25% (10) 


ez E,(E, — 28%) 


This form is appropriate among other things for studying kinematical properties 
of spin waves in regular antiferromagnetics superstructures (see e.g. Cofta 1960a). 
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